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Abstract 

Embedded random matrix ensembles are generic models for describing 
statistical properties of finite isolated interacting quantum many-particle 
systems. For the simplest spinless fermion (or boson) systems, with say m 
fermions (or bosons) in N single particle states and interacting via k -body in¬ 
teractions, we have EGUE(fc) [embedded GUE of k -body interactions) with 
GUE embedding and the embedding algebra is U(N). A finite quantum 
system, induced by a transition operator, makes transitions from its states 
to the states of the same system or to those of another system. Examples 
are electromagnetic transitions (then the initial and hnal systems are same), 
nuclear beta and double beta decay (then the initial and hnal systems are 
different), particle addition to or removal from a given system and so on. 
Towards developing a complete statistical theory for transition strength den¬ 
sities (transition strengths multiplied by the density of states at the initial 
and hnal energies), we have derived formulas for the lower order bivariate 
moments of the strength densities generated by a variety of transition oper¬ 
ators. Firstly, for a spinless fermion system, using EGUE(fc) representation 
for a Hamiltonian that is fc-body and an independent EGUE(t) representa¬ 
tion for a transition operator that is t-body and employing the embedding 
U(N ) algebra, finite-iV formulas for moments up to order four are derived, 
for the hrst time, for the transition strength densities. Secondly, formulas for 
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the moments up to order four are also derived for systems with two types of 
spinless fermions and a transition operator similar to beta decay and neutri¬ 
noless beta decay operators. In addition, moments formulas are also derived 
for a transition operator that removes ko number of particles from a system of 
m spinless fermions. In the dilute limit, these formulas are shown to reduce 
to those for the EGOE version derived using the asymptotic limit theory of 
Mon and French [Ann. Phys. (N.Y.) 95 (1975) 90]. Numerical results ob¬ 
tained using the exact formulas for two-body (k = 2) Hamiltonians (in some 
examples for k = 3 and 4) and the asymptotic formulas clearly establish 
that in general the smoothed (with respect to energy) form of the bivariate 
transition strength densities take bivariate Gaussian form for isolated finite 
quantum systems. Extensions of these results to bosonic systems and EGUE 
ensembles with further symmetries are discussed. 

Keywords: Finite many-particle quantum systems, Embedded ensembles, 
Transition strengths, Bivariate moments, U (. N ) Wigner-Racah algebra, 
Asymptotics, Bivariate Gaussian 
PACS: 05.30.-d, 21.60.Fw, 24.60.-k 


1. Introduction 

Wigner introduced random matrix theory (RMT) in physics in 1955 pri¬ 
marily to understand statistical properties of neutron resonances in heavy 
nuclei [1, 2], Depending on the global symmetry properties of the Hamil¬ 
tonian of a quantum system, namely rotational symmetry and time-reversal 
symmetry, we have Dyson’s tripartite classification of random matrices giv¬ 
ing the classical random matrix ensembles, the Gaussian orthogonal (GOE), 
unitary (GUE) and symplectic (GSE) ensembles. In the last three decades, 
RMT has found applications not only in all branches of quantum physics 
but also in many other disciplines such as econophysics, wireless communica¬ 
tion, information theory, multivariate statistics, number theory, neural and 
biological networks and so on [3-7]. However, in the context of isolated fi¬ 
nite many-particle quantum systems, classical random matrix ensembles are 
too unspecific to account for important features of the physical system at 
hand. One refinement which retains the basic stochastic approach but allows 
for such features consists in the use of embedded random matrix ensembles 

[8-13]. 

Finite quantum systems such as nuclei, atoms, quantum dots, small 
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metallic grains, interacting spin systems modeling quantum computing core 
and ultracold atoms, share one common property - their constituents pre¬ 
dominantly interact via two-particle interactions. Therefore, it is more ap¬ 
propriate to represent an isolated finite interacting quantum system, say with 
m particles (fermions or bosons) in N single particle (sp) states by random 
matrix models generated by random /c-body (note that k < m and most often 
we have k = 2) interactions and propagate the information in the interactions 
to many (m) particle spaces. Then we have random matrix ensembles in m- 
particle spaces - these ensembles are defined by representing the /c-particle 
Hamiltonian (H) by GOE/GUE/GSE and then the m particle H matrix is 
generated by the m-particle Hilbert space geometry. The key element here is 
the recognition that there is a Lie algebra that transports the information in 
the two-particle spaces to the many-particle spaces. As a GOE/GUE/GSE 
random matrix ensemble in two-particle spaces is embedded in the m-particle 
H matrix, these ensembles are more generically called embedded ensembles 
(EE). 

Embedded ensembles have proved to be rich in their content and wide in 
their scope. A book giving detailed discussion of the various properties and 
applications of a wide variety of embedded matrix ensembles is now available 
[13]. Significantly, the study of embedded random matrix ensembles is still 
developing. Partly this is due to the fact that deriving generic properties of 
these ensembles is not mathematically tractable and this is the topic of the 
present paper. A general formulation for deriving exact analytical results is 
to use the Wigner-Racah algebra of the embedding Lie algebra [13]. Till now, 
this has has been applied only in the study of one- and two-point functions in 
eigenvalues. The focus in the present paper is on transition strengths which 
on one hand probe the structure of the eigenfunctions of a quantum many- 
body system and on the other are important ingredients in many applications 
(for example, beta decay transition strengths are essential for nucleosynthesis 
studies). Also, as emphasized in [2, 8], there are many open questions in the 
random matrix theory for transition strengths in finite interacting quantum 
many-particle systems. 

For finite quantum many-particle systems, induced by a transition oper¬ 
ator, a given system makes transitions from its states to the states of the 
same system or to the states of another system. Examples are electromag¬ 
netic transitions (then the initial and final systems are same), nuclear beta 
and double beta decay (then the initial and final systems are different), parti¬ 
cle addition to or removal from a given system and so on. Given a transition 
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operator O acting on the m particle eigenstates | E) of H will give transi¬ 
tion matrix elements | (Ef \ O \ Ej) | 2 . Fig. 1 shows a schematic picture of 
transition strengths. In the statistical theories, it is more useful to deal with 
the corresponding transition strength density (this will take into account 
degeneracies in the eigenvalues) defined by 

Io(E„E f ) = I(E f ) \{E, | O | Si)| 2 I(E 0 . (1) 

In Eq. (1), 1(E) are state densities normalized to the dimension of the m- 
particle spaces. Note that E t and Ef belong to the same m-particle system 
or different systems depending on the nature of the transition operator O. 
In the discussion ahead, we will consider both situations. Random matrix 
theory has been used in the past to derive the form of the smoothed 1(E). 
In particular, exact (finite N) formulas for lower order moments (H p ), p = 2 
and 4 of 1(E) are derived both for EGOE and EGUE ensembles using group 
theoretical methods directly [14] or indirectly [15]. Let us add that some 
results valid only in the asymptotic limit (essentially N —> oo, m —> oo, 
m/N —> 0) are also available in literature for the state densities [16, 17]. 
Going beyond the eigenvalue densities, here we will apply group theoretical 
methods and derive for the first time some exact formulas for the ensemble 
averaged lower order moments of the transition strength densities Io(Ei, Ef). 
These then will give the general form of the smoothed (with respect to both 
Ei and Ef) transition strength densities. We will focus on fermionic systems 
and at the end, discuss the extension to bosonic systems. It is important 
to mention that in the present paper fluctuations in transition strengths are 
not considered and they will be addressed in a future publication. For some 
discussion on strength fluctuations see [2, 8, 9]. More importantly, using 
smoothed transition strengths statistical spectroscopy analysis of transition 
strengths (for example for particle transfer, electromagnetic transitions, beta 
decay and even double beta decay matrix elements) in nuclei and also in 
atoms is possible [18, 19]. In addition, smoothed form of the strength densi¬ 
ties can be used to calculate transition strengths that are needed for example 
in astrophysical applications [20, 21] and for neutrinoless double beta decay 
[19, 22, 23]. These formed the main motivation for the present study. Let us 
add that in the past, besides deriving the dilute limit formulas for the bivari¬ 
ate moments of the transition strength densities in some situations [24, 25], 
there are suggestions of using a polynomial expansion theory [26] (later in 
[24] it was shown that the polynomial expansion starts with the GOE result 
and hence in general inappropriate) for transition strengths, a specialized 
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theory for one-body transition operators [18, 27-29] and using the bivariate 
t-distribution form for transition strength densities [30]. 

In general, the Hamiltonian may have many symmetries with the fermions 
(or bosons) carrying additional degrees of freedom such as spin, orbital angu¬ 
lar momentum, isospin and so on. Also, the system may comprise of different 
types of fermions (or bosons); for example, in atomic nuclei we have protons 
and neutrons. In addition, a transition operator may preserve particle num¬ 
ber and other quantum numbers or it may change them. Of all these various 
situations, here we will consider three different systems in detail, (i) Firstly, 
we will consider a system of spinless fermions and a transition operator that 
preserves particle number, (ii) Secondly, we will consider a system with two 
types of spinless fermions with the transition operator changing /c 0 number 
of particles of one type to ko number of particles of other type as in nuclear 
beta decay and neutrinoless double beta decay, (iii) Third system considered 
is transition operators that remove (add) say ko number of particles from 
(to) a system of m spinlcss fermions. In all these we will restrict to EGUE. 
Before giving a short preview, let us mention that some of the results in the 
present paper are reported in three conferences earlier [31-33] and in a thesis 
[25]. . ' 

Section 2 gives some basic results for EGUE(/c) for spinless fermion sys¬ 
tems as derived in [14] and some of their extensions. Using these results, 
formulas for the lower order bivariate moments of the transition strength 
densities for the situation (i) above are derived and they are presented in de¬ 
tail in Section 3. The corresponding asymptotic limit formulas are presented 
in Section 4. Similarly, results for the situation (ii) above are given in detail 
in Section 5 and their asymptotics are given in Section 6. In addition, in 
Section 7, results for the situation (iii) are given along with the correspond¬ 
ing asymptotic formulas. Extensions to bosonic systems is discussed briefly 
in Section 8. Finally, Section 9 gives conclusions and future outlook. 

2. Basic EGUE(fc) results for a spinless fermion system 

Let us consider m spinless fermions in N degenerate sp states with the 
Hamiltonian H a /c-body operator, 

H = Y J V ij (k)A\(k)A j (k), V ii {k) = (k,i\H\k,j) . (2) 
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Here A\(k) is a k particle (normalized) creation operator and Ai(k) is the 
corresponding annihilation operator (a hermitian conjugate). Also, i and j 
are A;-particlc indices. Note that the k and m particle space dimensions are 
(^) and (] v ) respectively. We will consider H to be EGUE(A;) in m-particle 
spaces. Then V tJ form a GUE with V matrix being Hermitian. The real 
and imaginary parts of Vij are independent zero centered Gaussian random 
variables with variance satisfying, 

V ab (k)V cd (k) = V 2 5 ad 5 bc . (3) 

Here the ‘overline’ indicates ensemble average. From now on we will drop 
the hat over H and denote, when needed, H by H[k). Let us add that in 
physical systems, k = 2 is of great interest and in some systems such as 
atomic nuclei it is possible to have k = 3 and even k = 4 [34, 35]. 

The U(N ) algebra that generates the embedding, as shown in [14], gives 
formulas for the lower order moments of the one-point function, the eigen¬ 
value density 1(E) = ((S(H — E ))) and also for the two-point function in the 
eigenvalues. In particular, explicit formulas are given in [13, 14] for ( H p )' n , 
P = 2,4 and ( H p ) m (H Q ) m , P + Q = 2,4. Used here is the U(N ) tensorial 
decomposition of the H(k ) operator giving v — 0,1 ,.,, i: k irreducible parts 
B u,UJv (k) and then, 

k 

H(k)= W^(k)B^(k). (4) 

v=0;lJv£v 

With the GUE(/c) representation for the H(k) operator, the expansion co¬ 
efficients W’s will be independent zero centered Gaussian random variables 
with, by an extension of Eq. (3), 

W^Jk) = Vl . (5) 

For deriving formulas for the various moments, the first step is to apply the 
Wigner-Eckart theorem for the matrix elements of B l ' ,UJv (k). Given the m- 
fermion states | f m Vi), we have with respect to the U(N ) algebra, f m = {l m }, 
the antisymmetric irreducible representation (irrep) in Young tableaux nota¬ 
tion and Vi are additional labels. Note that v introduced above corresponds 
to the Young tableaux {2 u l N ~ 2u } and oj v are additional labels. Now, Wigner- 
Eckart theorem gives 

(fmVf | B^(k) | f mVi ) = (f m || B v {k) || f m ) . (6) 

JmVf , JmVz 
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Here, is the reduced matrix element and CZZ _is a U(N) 

Clebsch-Gordan (C-G) coefficient [note that we are not making a distinction 
between U(N ) and SU(N)]. Also, if | f m v t ) represents a m-particle state, 
then | f m vl } represents a m-hole state (see [14] for details). In Young tableaux 
notation f m = It is important to mention that \uuj u ) = \VH y). Fig. 

2 shows some typical Young tableaux and also the Young tableaux / that 
corresponds to a given {/}. Definition of B ,AUJv (k) and the U(N) Wigner- 
Racah algebra will give, 

I {fm. || B u (k) II f m ) | 2 = A v (N,m,m - k) , 


A M (Ad, m', r) 



^N' — m! + r — /j^j 


(7) 


Note that A u (N, m,k) is nothing but, apart from a N and m dependent 
factor, a U(N) Racah coefficient [14]. This and the various properties of 
the U (. N ) Wigner and Racah coefficients give two formulas for the ensemble 
average of a product of any two m-particle matrix elements of H, 


(fmV 1 | H(k) | f m V 2 ) (f m V 3 | H{k) | f m Vi) 


Vl V A V (N, m, m — k) C u / 

< ^ Jrr 


fm^l 5 fmV 2 fm^S ? fmV 4 5 


( 8 ) 


and also 


(fmV 1 I H(k) I f m V 2 ) (f m v 3 | H{k) | / m U 4 ) 

m—k 

fmV 1 ; fmV4 fmVs ? fmV 2 

I/= 0 \LJ u 


- V 2 

— V H 


T A “(N,m,k)CT 

< Jm 


(9) 


Eq. (9) follows by applying a Racah transform to the product of the two 
C-G coefficients appearing in Eq. (8). Let us mention some properties of the 
U(N) C-G coefficients that are quite useful in deriving the formulas given in 
Sections 3 and 5, 


• ^ Jm'Vi > Jm'Vi 

Vi 



C T T~ _ 

JmVi 5 JmVj 



& 


Vi ,Vj 
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(jfabVab _ /_1 )<t>(fa,fb,fab) (jfabV a b 

faVa , fbVb ^ Jb v b , faV a ' 


E s~lV,Uh, ,U> u l _ r r 

, JZvJ fmVj , JZvi VV ' 


( 10 ) 
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Here, (—1)^ is a phase factor depending on the irreps f a , ft, and f a b. See [14, 
36-38] for further details on the phase relations for U(N) C-G coefficients. 
From now on we will use the symbol f m only in the C-G coefficients, Racah 
coefficients and reduced matrix elements. However, for the matrix elements 
of an operator we will use m implying totally antisymmetric state for fermions 
(symmetric state for bosons). 

Starting with Eq. (4) and using Eqs. (5), (9) and (10) will immediately 
give the formula, 


(H 2 (k)) m 



El (mvi | H 2 {k ) | mvi) 

Vi 


Vl A 0 (IV, m, k) . 


( 11 ) 


Similarly, for (H 4 ) m , the ensemble average is decomposed into three terms 
as 


(H 4 (k)) m = E ( m A I H A (k) | mvi) 

Vi 

= (mvi | H(k) | mvj ) (mvj \ H(k) \ mv p ) (mv p \ H(k) | mvi) (mvi \ H(k) \ mvi) 


Vi,Vj,V p ,Vi 


= E 


(mvi | H(k) | mvj) (mvj \ H(k) \ mv p ) ( mv p \ H(k) \ mvi ) (mvi \ H(k) \ mvi) 


+ (mvi | H{k) | mvj) (mvi \ H{k) \ mvi) (mvj \ H(k ) | mv p ) (mv p \ H(k) \ mvi ) 


(mvi | H(k ) | mvj) (mv p \ H(k) \ mvi) ( mv j I H(k) \ mv p ) (mvi \ H(k ) | mn*)] . 

( 12 ) 


and the 


It is easy to see that the first two terms simplify to give 2 

third term is simplified by applying Eq. (8) to the first ensemble average and 
Eq. (9) to the second ensemble average. Then, the final result is 


9 /» r\ - 1 

{H\k)) m = 2 \{H 2 (k)) m ] +Vfi( ) A v (N,m,k)A v (N,m,m^k)d(N : v) , 

L J V m / iZo 

(13) 

where 


d{N : v) = 


2 / jy \ 2 

V ~ 1 


(14) 



Finally, a by-product of Eqs. (9) and (10) is 


y, ( mvi | H{k) | mvj) (mvj \ H(k ) | mv k ) = (H 2 {k)) m 5 VitVk , (15) 

v 3 

and we will use this in Section 3. Now we will discuss results for the moments 
of the transition strength densities generated by a transition operator O. 

3. Lower-order moments of transition strength densities: H EGUE(fc) 
and O an independent EGUE(t) 

For a spinless fermion system, similar to the k- body H operator, we will 
consider a t-body transition operator O represented by EGUE(f) in the m- 
particle spaces. Then, the matrix of O in i-particlc space will be a GUE 
with the matrix elements O a b(t) being zero centered independent Gaussian 
variables with the variance satisfying, 


O ab (t)O cd (t) = V$6 ad 6 bc . (16) 

Further, we will assume that the GUE representing H in /e-particle spaces 
and the GUE representing O in t particle spaces are independent (this is 
equivalent to the statement that O does not generate diagonal elements 
(Ei | O | Ej)\ see [24]). It is important to mention that in the past there 
were attempts to numerically study, in some nuclear (2sld) shell examples, 
transition strengths using the eigenvectors generated by a EGOE (with sym¬ 
metries) representation for H but taking O to be a realistic transition oper¬ 
ator. The results are found to be in variance with those obtained using a H 
defined by a realistic two-body interaction and O a realistic transition oper¬ 
ator [39]. As described in the present paper, a proper random matrix theory 
for transition strengths has to employ ensemble representation for both the 
Hamiltonian and the transition operator. 

With EGUE representation, O is Hermitian and hence, O t = O. Mo¬ 
ments of the transition strength densities Io(Ei,Ef) are defined by 

M PQ (m ) = (Crt(t)HQ(k)0(t)H p (k)) m = {0{t)H^(k)0(t)H p (k)) m . (17) 

Here the ensemble average is w.r.t. both EGUE(/c) and EGUE(f). Now we 
will derive formulas for Mpq with P + Q = 2 and 4; the moments with odd 
value of (P + Q) will vanish by definition. 
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Firstly, the unitary decomposition of 0(t) gives, 


o(t)= y. ■ ( 18 > 

iy=0;uj l/ eiy 

The [f’s satisfy a relation similar to Eq. (5). Now, for P = Q = 0, using Eq. 
(11), we have 

(Oit) 0{t)) m = Vq A°(N, m, t ) . (19) 

Moreover, we have the relations 

(0(t)0(t)H p {k)) m = ( 0(t)0(t)) m (H p (k)) m = (0{t)H p (k)0{t)) m (20) 

and their proof is as follows. Let us consider (0(t)0(t)H p (k)) m . Then, 

(■mvi | 0(t)0(t ) | mvj)(mvj \ H p (k ) | mvi) . 

Vi,Vj 

_ _ ( 21 ) 

Now applying Eq. (15) gives (mvi \ 0(t)0(t) \ mvj) = (0(t) 0(t)) m S VijVj . 
Substituting this into Eq. (21) will give the first equality in Eq. (20). The 
second equality (0(t)0{t)H p (k)) ru = ( 0(t)H p (k)0(t)) m follows from the 
cyclic invariance of the m-particle average. Equation (20) gives the moments, 


m)om F m m = ( A ) 

\m J 


M 20 (m) = M 02 (m) = (0(t) 0(t)) m (H\k)) m , 
M 40 M = M 04 (m) = (0(t) 0(t)) m (H\k)) m . 


( 22 ) 


Formulas for (O(t) O{t )) m , ( H 2 (k)) m and ( H 4 (k)) m follow from Eqs. (19), 
(11) and (13). Thus, the non-trivial moments Adpq for P + Q < 4 are M 41 , 
M 13 = M 31 and M 22 . 

ft is easy to recognize that the bivariate moment M n has same structure 
as the third term in Eq. (12) for ( H 4 (k)) m as the O and H ensembles are 
independent. Then, the formula for M n follows directly from the second 
term in Eq. (13). This gives, 

M u (m) = (0{t) H(k) 0(t ) H(k)) m 


V 2 V 2 

v O V H 


— l min(k,m—t) 

^ A u (N, m,t) A u (N, m, m — k) d(N : v) 
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This equation has the correct t k symmetry. Using Eqs. (23) and (11), 
we have formula for the bivariate correlation coefficient £, 


CM = 


Mu (to) 


min(t,m—k ) 

A U (N, m, k) A V (N, to, to. — t) d(N : v) 

v=0 


y/ M 20 (to) M 02 (to) 


A°(A^, to, t) A°(iV, to., fc) 


(24) 

Turning to M P q with P + Q = 4, the first trivial moment is M 13 = M 31 . For 
M 31 we have, 


M 31 (m) = ( 0{t)H{k)0(t)H\k )) 

^ TV' 


^ (mvi | d)(t) | my,-) (my,- | if(fc) | m»fc) (ran*, | O(t) \ mvi) (mvi \ H 3 (k) \ mvi) 


Vi,Vj,V k ,Vi 


N' 


y] {mvi | 0(t) | mvj) (mvk \ 0(t) \ mvi) {mvj \ H(k) \ mvk) {mvi \ H 3 (k) | mvi) . 


Vi,Vj,V k ,V t 


(25) 

In Eq. (25), the last equality follows from the fact that the EGUE’s repre¬ 
senting H and O are independent. The ensemble average of the product of 
two O matrix elements follows easily from Eq. (9) giving, 


fm v j 


( mi’i | 0(t) | mvj) (mv k \ 0(t ) | mvi) = Uj ^ A v (N,m,t) J2 C f^v i ,f^v i C f^ k , 

u =0 u>v 

(26) 

The ensemble average of the product of a H matrix element and II ‘‘ matrix 
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element appearing in Eq. (25) is, 


(■mvj | H(k) | mvk) {mvi \ H 3 (k ) | mvt) 


1 IK 

iQ V ■* 

(mvj \H(k)\ 

mv k ) (mvi | 

H(k) 

miy) (m«p | (A:) mu,) (mv q 

| f7(fc) | mvi) 

V 

V p ,V q 

(mvj | H{k) 

mv/-) ( mvi 

I 7/(7-) 

| ?nr p ) (miip | H(k) mv q ) (mv g 

| 7f(fc) | mvi) 


+ (mvj | H(k) | mvk) (mv q \ H(k ) | mv^ (mvi \ H(k) \ mv p ) (mv p \ H(k ) | mv q ) 


+ (maj | i7(/c) | mrfc) (mvp \ H(k) \ mv q ) {mvi \ H(k) \ mv p ) (mv q \ H(k) \ mvi )J . 

(27) 

The first two terms in Eq. (27) simplify to give 2 (H 2 (k)) m Mn(m), using 
Eq. (26). Simplifying the third term using Eqs. (8) and (9), we get 


(mvj | H(k) | mvk) (mv p \ H(k) \ mv q ) (mvi \ H(k ) | mv p ) (mv q \ H(k ) | mvi) 


k m—k 


i/i=0;a; Vl i/2=0;cu^ 2 


(28) 

Combining this with Eq. (26) and applying the orthonormal properties of 
the C-G coefficients will give the final formula for M 31 , 


M 31 (m) = (0(t) H(k ) 0(t) H 3 (k)) m 


= 2 {{H(k)] 2 r M u (m) + 





min(k,m—k,m—t) 

A V (N, m, t) A V (N, m, k) tv v {N, m, m 

z/=0 


k) d(N : v) . 


(29) 


Thus, M 31 (m) involves only the A functions. 

To derive the formula for M 2 2 , we will make use of the decompositions 
similar to those in Eqs. (25) and (27). Then it is easy to see that M 22 will 
have three terms and let us say they are Si, S 2 and S 3 . The first term is 
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S 1 = (0(t)0(t)) m (H 2 (k)) 


l2 


and the second term So is 


N 


-l 


5 2 = r 1 S 2 i(v i ,v j ,v r ,vi) S 2 2(vi,v j ,v l ,vp,v Q ) ; 

,771/ z ' 

Vi,Vj,V r ,Vl,Vp,VQ 


S 2 i(vi,vj,v r ,vi ) = (m,Uj | £>(t) | m,u,) (m,ry | £>(£) | m,uj) , 


(30) 


S 22 (v i ,v j ,vi,v P ,v Q ) = (m, Vj | H(k) \ m,v P ) ( m,v Q \ H(k ) | m,Vi) 


x {m,vp | //(£;) | m,v r ) {m,vi \ H(k) \ m,VQ ) . 

The sum involving S 22 is simplified by applying Eq. (9) to the two ensemble 
averages in S 22 and using the orthonormal properties of the C.G coefficients. 
Now, applying Eq. (8) to the ensemble average in S 2 \ and then multiply¬ 
ing with S 22 will (after using again the orthonormal properties of the C-G 
coefficients) lead to the final result for S 2 , 

_1 min(t,m—k) 

J2 a V (N, m, m - t ) [A V (N, m, k)f d(N : u) . (31) 

u=0 

Similarly, the term S :i is 

^ = ( J S S 3l( v i> V V l, r , V h V P , Vq) S 32 (Vi, V r , V P , Vq) ; 

' ' Vi,Vj,V r ,Vl,Vp,VQ 


S3i{v i ,v j ,v r ,vi,v P ,v Q ) = ( m,Vi I 0(t) \ m,Vj ) (m,v r \ 0(t ) | m,vi) 



x (m, Vj | H(k) | m,v P ) ( m,v t \ H(k) \ m,v Q ) , 


S 32 (vi,v r ,v P ,v Q ) = (m, v P | H(k) \ m,v r ) (m,v Q \ H(k) \ m,Vi ) . 

(32) 

The term S 3 1 is simplified by first using Eqs. (4) and (18) and this gives, 
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after carrying out the ensemble averages, 


t k 

S>. = VJV2 Y, EE 

Vi,Vj,Vr,Vl,Vp,VQ V\=0 ]L0 V1 !/ 2 =0;O )v 2 

X ( m,v r | B' ,lUu i{t) | to,i>/) (to, y,- | B V2 ^{k) \ m,v P ) ( m,vi \ B V2U ^(k) \ m,VQ) 

t k 

= ^ E E E (m, Vi | B Vl ^{t) B V2 ^{k) | m,u P ) 

Vi,V r ,Vp,VQ Vl=0',U) l , 1 V2=0\UJ„ 2 

x (m,v r | B uiu} "i (t) B V2U ^ (k) | m,VQ) . 

(33) 

Now, coupling the H’s in each matrix element in the last equality in Eq. (33) 
applying the Wigner-Eckart theorem for the matrix elements of the coupled 
tensor operator [B Vl {t)B U 2 {k)] UUv along with the application of Eq. (9) to 
the 632 term will give the formula for the S 3 term. With this, the M 2 2 is 
given by, 

m 22 { m ) = (o{t)o(t)) m [wwmr ] 2 

( v _i [ min(t,m—k) 

^ A u (N, m, m-t) [A u (N, m,k )} 2 d(N : u) 
m ' y u=o 

min(k-\-t,m—k ) t k 

+ Y, A v {N,m,k)d{N:v) YY J Y,\^\\[ BVi m u 2 m V ' P \\f m )\ 2 

v =0 ^i=0 i^ 2—0 p 

(34) 

Here p labels multiplicity of the irrep v in the Kronecker product v\ x v 2 —> v. 
We will see in Section 4 that Eq. (34) is useful in deriving asymptotic results. 
In order to evaluate M 2 2 , we need a formula for the reduced matrix element in 
Eq. (34). This is obtained by considering the corresponding matrix element, 
decoupling the coupled operator, using complete set of states between the 
two operators, applying Wigner-Eckart theorem to the two matrix elements 
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and then simplifying the sums over all the C-G coefficients. This gives, 


tf - lr(W 'l (35 

x (f m II || f m ) (f m || B^(k) || fm) f m v) p . 

Eq. (35) gives the formula for M 22 in terms of U (.AT) Racah coefficients, 

M 22 (m) = (0(t) H 2 {k) 0(t) H 2 {k)) m = V%V& { [A 0 (TV, m, k)} 2 A 0 (TV, m, t) 


+ 



-1 


min(t,m—k) 

^ A v (N,m,m — t) [A" (TV, m, &)] 2 < 7 (TV : u) 

u =0 


+ 



-2 


min(k-\-t,m—k) t k 

E EE A" (TV, m, /c) A*' 1 (TV, m, m — t) 

p =0 ui =0 1/2=0 


x Ar(N,m,m - k)d(N : !x 1 )d(A r : v 2 )^^[U(f m v 1 f m v 2 ; /m")J 2 

(36) 

The U or Racah coefficient in Eq. (36) is with respect to C/ (TV) and a formula 
for this is not available in closed form in the literature. Deriving formulas for 
this [/-coefficient is an important open problem. We will show ahead that it 
is possible to derive a formula valid in the asymptotic (TV — > 00, m —> 00, 
m/TV —* 0 and k,t fixed) limit. 



4. Asymptotic results for the bivariate moments for H EGUE(fc) 
and O an independent EGUE(f) 


Lowest order (sufficient for most purposes) shape parameters of the bi¬ 
variate strength density are the bivariate reduced cumulants of order four, 
i.e. k rs {m ), r + s — 4. The k rs {m ) can be written in terms of the normalized 
central moments Mpg(m) where Mpg(m) = MpQ(m)/M 00 (m). Then, the 
scaled moments /ipg(m) are 


/W m ) 


M pq (m) 


"...- 

P/2 

- -- 

Q/2 

M 20 (m) 


M 02 (m) 



-P + Q — 2 • 


(37) 
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Note that /r 2 o (m) = )Mn(jn) = 1 and p u {m) = £(m). Now the fourth order 
cumulants are, 


k 40 (m) = Haq (m) - 3 , k 04 (m) = /i 04 (m) - 3 , 
k 31 (m) = n 31 (m) -3£(m), k 13 (m ) = AM 3 (m) - 3£(m), (38) 

^22 (nr) = /r 22 (m) - 2 £ 2 (m) - 1 . 


The |/c rs (m)| ~ 0 for r + s > 3 implies that the transition strength density is 
close to a bivariate Gaussian (note that in our EGUE applications, k rs (m ) = 
0 for r + s odd by definition). Numerical results for k rs (m), r + s = 4 and also 
for £(m) for some typical values of N, m, k and t are shown in Table 1. It is 
seen that for sufficiently large values of N and m and k +1 relatively small, 
in general the magnitude of the fourth order cumulants is very small (< 0.3) 
implying that for EGUE, transition strength densities approach a bivariate 
Gaussian. However, with increasing k + t value it is seen that £(m) —» 0 and 
this is the GOE result. Also, in this limit the marginal densities approach 
semi-circle form giving k 4 o(rn) —> —1. For a better understanding of these 
results, it is useful to derive expressions for npQ(m) and thereby for kpQ^rn), 
using Eq. (38), in the asymptotic (asyrnp) limit defined by N —> oo, m —>■ oo, 
rn/N —> 0 and fixed k and t with t < k. First we will consider N oo and 
m fixed with k,t « m. Two relations we use are, 


f N — p\ p /N —¥ 0 N r 
\ r J r! 


d(N : v] 


v/N-tO N' 


2 v 


(u\) 


M 2 


(39) 


Note that d,(N : u) is defined in Eq. (14). Let us start with the formula for 
£ given by Eq. (24). Applying Eq. (39), it is easily seen that only the term 
with v — t in the sum in Eq. (24) will contribute in the asymptotic limit. 
Using this and applying Eq. (39) to the formula for A u given by Eq. (7) will 
lead to, 


£(m) 


m\N» (V)( 


m—t\ (N—m+k—t\ (N—2t 


HIT) 


(t!) 2 (™^ (N —' 1 ^ N—m+t ^ 


asyrnp (m 
k 


-i 


Similarly, /i 40 (m) will be, using Eqs. (22), (13) and (11), 
AUo (m) - 


o m\ N 2k ( m k k ) ( N k m ) ) asyrnp o 

T»T_ / 7 i \ o / \ /AT ... i 7,\ / ... \ /AT ... i I,\ ' " I 


JV” («) 2 afTUDfTU 


m — t 
k 


(40) 


-1 


m — k 
k 

(41) 
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Turning to /r 31 (m), it is easy to see from Eq. (29) that M 31 (m) has two 
terms. The first term is 2 £(m) (0(t) 0(t)) m [{H{k)H(k)) m } 2 and in the sum 
in the second term, only v = k will survive in the asymptotic limit. These 
then will give, 


Ahn (m) 


2 £( 777 ) 



1 A k (N, to, t)A k (N , m, k)A k (N , m, m — k) d(N : k ) 

A 0 (AT, 777, t)[A°(A[ 777, fc)] 2 


asymp 




2 + 


m \ 1 /m — k 


= £(m) fJ,4o{m) . 


(42) 

Here, Eq. (39) is used in the final simplihcations. Finally, let us consider 
h- 22 (ui). Firstly, M 2 2 (m) has three terms as seen from Eqs. (36) and (34) and 
let us call them XI, X2 and X3. Then, there will be corresponding three 
terms in p 22 {m) and we call them Tl, T 2 and T 3. It is seen that T1 = 1 
and T 2 is (in the corresponding X2 sum, only is — t term will contribute in 
the asymptotic limit), 


T 2 



A 4 (AT, 777 , ?77. — f) [A 4 (X, ?77, fc)] 2 d(X : t ) 
A 0 (AT, m, t)[A°(X, 777 , fc)] 2 




(43) 

As a formula for the [/-coefficient appearing in Eq. (36) is not available, we 
use Eq. (34) for simplifying T 3. Then, in the X 3 sum only is = t + k term 
will survive in the asymptotic limit giving 


(N\ - 1 A 4+fc (X, 7n,A:)d(X:t + fc) 
\777/ A 0 (AT, 777, t)[A°(N, 777, A))] 2 


(/ m || [B‘(t)i3‘(7)] ,+ ‘ 


For further simplification of T3, we will use the relation 



PQ {0{t)H{k)0{t)H{k)) m = Vq V »EEE I (fm 

' ' ui—0 U 2 —O u;p 

Then, T 3 in the asymptotic limit will be 


[B*(t)B v *(k)r P \\fm )\ 2 d(N : v) . 

(45) 


A t+k (N, 777 , k) 
A °(N,m,k) 


[C(X 00 )] 


fm\ 2 fm — t — k\ fm — t\ 

\k) V k ){ k )■ 

(46) 
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Notice that, here we have used Eq. (40) for £(m) in the N —> oo limit. Now 
combining Ti, T 2 and T 3 we have 


asymp 

/x 2 2(nr) —> 1 + 


—2 / ,\ 2 
m — V 


+ 


-2 


m — t — k 
k 


m — t 


(47) 


Comparing Eq. (46) with the formula for T 3 as given by Eq. (36), it is easy 
to see that in the asymptotic limit 


[u(.fm t fmk; fm. t + k)} 


2 asymp 


-i 


m — t 
k 


(48) 


The asymptotic formulas derived from the exact results and given by Eqs. 
(40), (41), (42) and (47) are same as those given before in [13, 24] where 
the theory given by Mon and French [16], that gives directly the asymptotic 
results, was used. This is a good test of the derivations in Section 3. It is 
important to add that in the derivations given in Section 3 we have considered 
only the binary correlated terms in the various sums [the binary correlations 
come in the ensemble average of the IE’s and C/’s defined by Eqs. (4) and 
(18)]. Within this constraint, the results in Section 3 are exact. On the 
other hand, in the Mon and French theory, the binary correlated terms are 
evaluated using formulas that are valid only in the asymptotic limit. 

Asymptotic formulas for the bivariate moments will give the asymptotic 
formulas for the fourth order bivariate cumulants /c 40 (m), k 3 i(m ) and k 2 2 (nr) 
by applying Eqs. (37) and (38). Then we have, 


ho(m) 

= k 0 4(m) 

fm — k 

\ fm 

= ( * 

)U 

k 31 (m) 

= k 13 (m ) 

= f(m)k 40 (m), 


fe 2 (nr) 

| 

fm - 


= £ 2 M < 

l( 


1 



Remember that £(m) is given by Eq. (40). Now, the 1/m expansion of k rs (m ) 
shows that k rs (m) = — k 2 /m + 0(l/m 2 ). This clearly demonstrates that we 
have for the transition strength densities bivariate Gaussian form in general 
in the dilute limit defined byiV—> 00 , m — > 00 , m/N —>■ 0 with fixed k and 
t. However, £(m) —> 1 in the dilute limit and this gives a singular bivariate 
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Gaussian that is unphysical. Therefore, in practice the dilute limit condition 
will not be realized and there will be departures from the bivariate Gaussian 
form. Thus, the corrections due to k rs (m), r + s — 4 should be included 
and out of the many ways to add the corrections, the Edgeworth expansion 
is considered to be the best [19, 40, 41]. The Edgeworth corrected bivariate 
Gaussian form including k rs (m ) with r + s — 4 is given in Appendix A. 

5. EGUE results for moments of transition strength densities: two 
types of spinless fermions with beta and double beta decay type 
transition operators 

Here we will consider a system with two types of spinless fermions with 
mi fermions of type #1 in Ni sp states and m 2 fermions of type ^2 in N 2 
sp states with the system H operator preserving (mi, m 2 ). This is similar to 
protons (p) and neutrons (n) in a atomic nucleus. Thus, we have two orbits 
with the first one having N\ sp states and the second N 2 sp states. Then the 
H operator, assumed to be £>body, is given by, 

H(k)= Ka:/3b(hj) A i(i) A p(i) A\(j) A b (j) ; 

i+j=k a, fiei a,bej 

Vaa:pb(i,j) = (l, Oi \ j, (1 \ H \ i, (3 \ j, 6) . 

Here, we are using Greek labels a, (3,... to denote the many particle states 
generated by fermions occupying the first orbit and the Roman labels a,b,... 
for the many particle states generated by the fermions occupying the second 
orbit. Note that A^i) creates the state | i,a) with i number of fermions in 
the first orbit and A' a (j ) creates the state | j,a) with j number of fermions 
in the second orbit. Similar is the action of the annihilation operators Ap(i) 
and A b (j). For example for a two-body Hamiltonian, (i,j) = (2,0), (1,1) 
and (0, 2). This corresponds to H = H pp + H pn + H nn for atomic nuclei; with 
‘p’ denoting protons and hi' denoting neutrons. For each (i,j) pair with 
i + -j — k, we have a matrix V(i,j) in the /c-particle space and H matrix 
is a direct sum of these matrices in the k particle space. Their dimensions 
being ( , ^ 1 ), N X N 2 and (^ 2 ) respectively. Action of the H operator on the 
\mi,v a : m 2 ,v a ) states of a (mi,m 2 ) system generates (mi,m 2 ) particle H 
matrix; v a and v a are respective additional labels. The H matrix dimen¬ 
sion is d(mi,m 2 ) = (^) (^ 2 ). To proceed further, the V(i,j) matrices are 
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represented by independent GUEs with matrix elements being zero centered 
Gaussian variables with variance, 


baa:/3b(b j') ba , a , :/3 , fe / (^ ij ) bj^(i, j) Jjj' ^jj' ^a/3' ^ab’ 3a’ &ba' ■ (51) 

It is important to note that the embedding algebra for the EGUE generated 
by the action of the V ensemble on \mi,v a : rn 2 , v a ) states is the direct sum 
algebra U(N i) © U(N 2 ). Thus we have EG\JE(k)-[U (Nj © U(N 2 )] ensemble. 

Going beyond the simple operators considered in Sections 2 and 3, here 
we will consider the following transition operator, 

O(ko) = ^2 O aa Al(k 0 ) A a (k 0 ) ; O aa = (k 0 , a \ O \ k 0 , a) . (52) 

a,a 

It is easy to see from Eq. (52) that the operator O(k 0 ) changes k 0 number of 
fermions in the second orbit to A; 0 number of fermions in the first orbit. This 
is similar to the action of /3 decay operator (then ko = 1) and neutrinoless 
double beta decay operator (/c 0 = 2). Eq. (52) gives the matrix elements 
of O for the matrix representation of O in the defining space and this is in 
general a rectangular T>a x V B matrix where V A = (^) and V B = (^). 
Just as in Section 3, we will assume GUE representation for the O matrix in 
the defining space and then, 

Oaa 0/31 = V£ S a/ 3 5 a b . (53) 

Note that the O matrix in the defining space is a GUE in the sense that the 
real and imaginary parts of the O aa are zero centered independent Gaussian 
variables with variance given by Eq. (53). With the GUE representation for 
the O matrix in the defining space, we have EGUE for O in many particle 
spaces. This matrix will be again a rectangular matrix with matrix elements 
connecting (mi,m 2 ) states with (mi + ko, m 2 — ko) states by the O operator. 
Fig . 3 shows an example for the H and O matrices in the defining and 
many-particle spaces. With independent GUE representations for H and O 
matrices in the defining spaces (then V aa :/ 3 b(i,j) are independent of O a ’ a '), 
we will derive formulas for the bivariate moments of the transition strength 
density 

I^ m2 \E h E f ) = /(m 1 +fco.m 2 -fco)( jB/ ) |<( mi + k 0 ,m 2 - ko),E f I o I (murmlEi)| 2 . 

(54) 


20 



Then, the ensemble averaged bivariate moments are 


M PQ (m 1 ,m 2 ) = (0^k 0 )HQ(k)0(ko)Hr(k)) {mi ’ m2) . (55) 

Note that O takes (mi, m 2 ) to (mj + /c 0 , m 2 — /c 0 ) uniquely and therefore the 
later is not specified explicitly in Eqs. (55) and (54). Also, it is important 
to note that O t ^ O. 

Given one type of spinless fermions, we easily have 



Now, substituting complete set of state between the and A operators and 
applying the Wigner-Eckart theorem [note that A\{k) and Aj(k) transform 
as the U(N) tensors f k = {l fc } and {f k } = {1^^°}] will give, 


(m || A^(k 0 ) || m- k 0 ) (m - k 0 || A(k 0 ) || m) 


G ' 1 °) • 

\m - k 0 J 


(m || A(k 0 ) || m + k 0 ) (m + k 0 || A\k 0 ) 



(57) 


For the derivations given ahead it is important to recognize that A^ (ko) and 
Ap(k 0 ) transform as the U(Ni) tensors f ko = {l fc °} and {f ko } = {l^ 1- ^’ 0 } 
and similarly A' a (ko) and A^iko) with respect to {7(1V 2 ). Moreover, using Eq. 
(56), we have 


o(k 0 )) (m '™> = vg “ mi ) , 

(58) 

(o(k 0 )ot(k 0 )) im ‘™> = vg ( iV 2 ( r m2 ) ■ 

Before turning to the bivariate moments MpQ(mi,m 2 ), let us consider 
the second and the fourth moment of the state densities I mi,rn2 (E). For 
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deriving formulas for these moments, unitary decomposition of H(k ) needs 
to be carried out with respect to U(N 1 ) © U(N 2 ) algebra and this gives the 
B v,Wv (k') tensors in U(N\) and U(N 2 ) spaces. In order to distinguish the 
tensors, we call them B and C in the two spaces respectively. Denoting 
these tensorial ranks by U\ and v 2 respectively, we have 

i+j=k i^i=0 i/2=0 

(59) 

Just as in Section 2, here also it can be proved that the W s will be indepen¬ 
dent zero centered Gaussian variable with variances satisfying, 

Wp’ (y \, Wj/j , V'2, 0J U2 ) l l r / (V| . U© ) ©, C©' ) 1// &i, r &j,l Gv, -V/ ‘ 

( 60 ) 

Using Eqs. (59) and (60) and the results given in Section 2, it is straight 

forward to derive the formulas for (H p (k))^ rni,rn2 ^ for P = 2 and 4. We have 
for P = 2, 

= V V|(i,j)A 0 (JV 1 ,m 1 ,i)A l, (JV 2 ,m 2 ,j). (61) 

Similarly for P = 4, 

<^(fc)> (mi ’ m2) = 2 [ (i/ 2 (U) (mi ’ m2) ] 2 + E E Va(i,3)V^(i',j') X(N 1 ,m 1 ,i,i')Y(N 2 ,m 2 ,j,j')-, 

i+j=k i'+j'=k 

( , _1 min(i,mi-i') 

M ^ A i/1 (ATi, toi, mi — *) A vl (Ni, toi, *') d(W : Vi) , 

1 / 1=0 

Y(N 2 , m 2 ,j, f) = ( 2 ) V A 122 (AT 2 , m 2 , m 2 - j) A 122 (7V 2 , m 2 ,j') d(N 2 : v 2 ) . 
\j n 2j ^ 

(62) 

It is also easy to show using Eqs. (56) and (53) that, 
(O^k 0 )O{k 0 )H p (k)) {mum2) = (O^k 0 )O(k 0 )) imi ’ m2) (H p (k)) {mi ’ m2) , 

{O^k 0 )H p (k)O(k 0 )) imi ’ m2) = {O^k 0 )O(k 0 )) {mi ’ m2) (HP(k)) {mi+ko ’ m2 - ko) . 

(63) 
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The bivariate moments M 00 , M 2 o, M 02 , M 40 and M 04 follow directly by ap¬ 
propriately using Eqs. (58), (61)-(63) in Eq. (55). 

Clearly, the first nontrivial bivariate moment is Mu (mi, m 2 ) and it is 
given by, 


Mu(mi,m 2 ) = (O^k 0 )H(k)O(k 0 )H(k)) m = ( ) 

L \ m iJ \rn2J J 

X ^2 Eavg{(m 1 ,a 1 ;m 2 ,a 1 | 0\k Q ) | m 1 k 0 ,a 2 ;m 2 - k 0 ,a 2 ) 

ai,a2,ct3, a 4'-ai,a2,CL3,a4 


x (m-i + k 0 , a 2 ] rn 2 - ko, a 2 \ H(k ) | mi + k 0 , a 3 ; m 2 - &o, ^ 3 ) 
x (m 4 + fe 0 , a 3 ; m 2 - fc 0 , a 3 | O(k 0 ) | m 4 , cn 4 ; m 2 , a 4 ) 
x (mi,a 4 ; m 2 , a 4 | //(/c) | m 4 , a 4 ; m 2 , a 4 )} 
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X ^2 Eavg { (r77-i, ap m 2 , % | O\k 0 ) | m x + k 0 , a - 2 ; m 2 - k 0 , a 2 ) 

cx. i , 0:2 ,c*3 ,a4 :ai ,d 2 ,a3 ,a4 


x (mi + k 0 , a 3 ; m 2 - ho, a 3 \ O(k 0 ) | mi, a 4 ; m 2 , a 4 )} 
x Eavg {(mj + fc 0 , <*2; - fe 0 , a 2 | H(k) \ mi + k 0 , a 3 ; m 2 - k 0 , a 3 ) 


x (mi,a 4 ;m2,a 4 | if(fc) | mi, ay m 2 , a 4 )} 


V 2 
v o 


(N ! 
V m i 



}"E E 

^ ai,a2,Q!3,Q!4:cti,a2,a3,c!.4:Q:,a 


U,OJ U) v' ,bl v t 


x (m 4 , ai; m 2 , a 4 | A* (k 0 )A a (k 0 ) \ m x + fc 0 , a 2 ; m 2 - fc 0 , a 2 ) 
x (mi + k 0 , a 3 ; m 2 - k 0 , a 3 | Ht(£;o) At(M I mi, a 4 ; m 2 , a 4 ) 
x (mi + k 0 ,a 2 ,m 2 - k 0 ,a 2 \ B v ' Uv (i)C u '’ w "’ (j) \ m i + k 0 ,a 3 ,m 2 - &o,a 3 ) 
x (mi,a 4 ;m 2 ,a 4 | B v,Uv {i)C v '' Wv ' (j) | mi, oq; m 2 , ai) . 

(64) 

Here we have used ‘Eavg’ to denote ensemble average instead of using 
the ‘overline’. We will use this notation at a few other places in the re¬ 
minder of this paper. In the second step above we have introduced com¬ 
plete set of states between the operators to write the ensemble average of 
((W (k 0 )H(k)O(k 0 )H(k)) as the ensemble average of four matrix elements. In 
the next step, this ensemble average is written as the product of the ensemble 
average of O matrix elements and the ensemble average of H matrix elements 
using the independence of the O and H ensembles. Finally we have used Eq. 
(59) to get the final form in Eq. (64). Now, it is easy to see that M u (rrii, m 2 ) 
factorizes into product of terms in mi and m 2 spaces giving, Mn(mi,m 2 ) to 
be of the form Vg [(^) ©] _1 ^ 2 i+j=k V£(i,j)X( N u m u i, k 0 )Y(N 2 , m 2 , j, k 0 ) 
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where for example X(Ni,rrii,i,k 0 ) is 
X(Ni,mi,i, ko) = 


y: {m i,ai 

a .1 


A a {ko) | mi + ko, a 2 ) (m\ + &o, a 3 | ^(^0) I ' m i> 04) 


x (m 1 + ko, a 2 \ B v,,J}v {i) \ rrii + k 0 , a 3 ) (mi, a 4 | B v,Wv (i) \ mi, ai) . 

. (65) 

Similarly Y(N 2 ,m 2 ,j,ko) can be written. Now, applying the Wigner-Eckart 
theorem along with Eqs. (7) and (57), we will be left with four C-G coef¬ 
ficients. These can be simplified to give a U(Ni) U— or Racah coefficient. 
Similarly, for the Y function we will get a U(N 2 ) Racah coefficient. Putting 
these together will give the final formula for Mu and for later usage we will 
write it in the following form, 


Mn(mi,m 2 ) = Vg 

i 

Y Xu(Ni,mi, k 0 ,i, v{) 


NX fN 2 
mi) \m 2 


y v H(hj) 


i+j=k 


Ni-ko 


x 


.1/1=0 


y Yu(N 2 , m 2 , k 0 ,j, u 2 ) 
R 2=0 


(N 2 - ko\ 
\m 2 - k 0 ) 


Xn(Ni,mi,k Q ,i,u) = 


l :) ■ 


1/2 


d(Ni : u) 

x [A u {N u m u mi — i ) A u (Ni, mi + ko, mi + ko — i )]^ 2 

X (_l)Wm 1 + k 0 ,/» 1 ,/* 0 )+W» 1 ,/» 1 ,>') U(f mi+ko y i f mi+ko fm 1 ; fko v) , 


Yu{N 2 , m 2 , ko,j, v) = Xu(N 2 , m 2 - k 0 , k 0 ,j, v) 



1/2 


d{N 2 : v) 


[A v {N 2 , m 2 , m 2 - j ) A U (N 2 , m 2 - k 0 , m 2 - k 0 - j )} 1/2 


x (_l)^(/ m2 -/ m2 -fcoTfco)+^/-2-fco-/ m2 -fc 0 .G U(f m2 / m2 _fc 0 f m2 fm 2 -k 0 ; fko v) . 

( 66 ) 

Note that f r = {l r } and u = {2 U , \ N ~ 2v } with N = Ni 01 N 2 as appropriate. 
Also, the [/-coefficient in Xu is with respect to U(Ni) while the one in Yu is 
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with respect to U(N 2 ). Formula for the [/-coefficient appearing in Eq. (66) 
is available in [36] and it is given by 

U({l n } {l n } {P} ; {l n_p } {2 u 1 n ~ 2u }) = (—1 )<H JV > n >P> 1 ') 

(p!) 2 (n — z/)! (AT + 1)! [(N — n)\] 2 (N — v — p) \ (N — 2is + 1) 1 1/2 
(z/!) 2 (p — u)\ (n — p)\ [(At + 1 — z/)!] 2 (AT — n + p)\ (AT — n — v)\ 

(67) 

The phase factor (—1)^ in Eq. (67) depends on the phase convention [36] 
and we will fix this phase later. Eqs. (64)-(66) show that in evaluating 
the bivariate moment (O^(k 0 )H^(k)O(k 0 )H p (k)) mi ' m2 we can use H(k) as 
Yli+j=kHi(i)H 2 (j) and then the moment will be a sum of terms where each 
term is a product of two functions with one in the rri i space [generated by 
Hi(i) with body rank i] and other in the m 2 space [generated by H 2 (j) with 
body rank j], Also, these functions follow from the results in Sections 2 and 
3 by appropriate application. This will be seen in the formulas for the fourth 
order moments that are discussed below. 

Turning to the fourth order moments, we need M 13 , M 31 and M 22 (we 
have already discussed M 40 and M 04 ). As ^ O , here M 13 ^ M 31 and 
similarly M 40 ^ M 04 (also M 20 ^ M 02 ). Following the procedure used for 
deriving the formula for Mu(rrii,m 2 ) and the results [Eqs. (25) and (27)] 
given in Section 3 for M 34 (m), we have for M 3 i(mi,m 2 ), 

M 31 (m 1: m 2 ) = ( Oi(k 0 )H(k)O(k 0 )H%k)) mi ’ m2 

= 2 Mn(m u m 2 ) + V’ | (^) } 

- E E V H (' h , ji) V H (A, j2) X 31 (Ah, m 4 , A, i 2 , k 0 ) Y 31 ( N 2 ,m 2 , ji , J2, k 0 ) ; 

h+ji=k h+j2=k 
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X 3 i(N 4 , mi, *i, * 2 , /c 0 ) = Y (mi, a t \ A a (k 0 ) \ mi + k 0 , a 2 ) 

ai,a2,(ys,a4,a5,aQ:a,i'i,uji> 1 ,V2 ->i*>v 2 

x (mi + k 0 , a 3 | AJ,(fc 0 ) | m u a 4 ) (m i + fc 0 , a 2 | B ui ’“ v i(i 1 ) \ mi + k 0 , a 3 ) 
x (mi,a 5 | B ui ’ u ^(ii) \ m 1 ,a 6 ) (mi ,« 4 | B V 2 ^(i 2 ) \ m 1 ,a 5 ) 
x (mi,a 6 | B l " 2 ' Uu ‘z (i 2 ) | my ay) , 


*31 (W 2 , m 2 ,ji,j 2 , fa) = Y ( m 2, bi | At (fco) | m 2 - fc 0 ,6 2 > 

bl,fe,fe3,f>4,65,b6:a,l'3A' ! ,3,I'4,Wi,4 

x (m 2 - A: 0 , &3 | A a (/co) I m 2 , & 4 ) (m 2 - & 0 , b- 2 | I m 2 - k 0 , b 3 ) 

x (m 2 , 6 5 | C U3 ’ UJl '3(j 1 ) | m 2 ,& 6 ) (m 2 ,b 4 \ C v ^ v *(J 2 ) \ m 2 ,b 5 ) 


x (m 2 , b 6 | C U4 ’ UJ ^(j 2 ) \ m 2 ,b 1 ) . 

( 68 ) 

The term X 3i (similarly Yu) is simplified using Eqs. (6), (7), (9) and (57) 
giving, 


X 3 i(Ni,mi,ii,i2, k 0 ) = 


Ni-ko 

mi 


i\ mi —22 

E E 

v\ =0 122=0 


A^TVymyi,) 


x [A^(A^i,mi, mi - i\) A vl (Ni,mi + k 0 ,mi + k 0 - fi)] 1/2 


x 


E 


C 


cy 


k'l 


V 2 ,Uv 2 


V\ ,uj v1 

fm\+k§° l -2 5 /mi+fcQ a 3 /m^5 5 /m^6 /m^4 1 /mi 0^1 


C, ’ 


Oi 1 ,Q!2,Q;3,Q!4,^5,CK6-CK 


(69) 


x C. 2 ^-_C 


(jko,OL 


fm-^OLQ , fm-^Ot .5 , /mi +/cq <^2 /mj+fcQ a 3 > fm^ OL 4 

The sum over a 5 and a 6 of C?’“^ C^Z 7— wil1 S ive <W 2 <U lA ,„ 2 - 

Jmi 0:5 , Q!6 jm^e i Jm^S 1 z 

Now, the remaining four C-G coefficients sum up to give a Racah coefficient. 
Carrying out a similar simplification of the C-G coefficients in Y 3i , we finally 
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obtain, 


M 3 i(m u m 2 ) = {Of(k 0 )H(h)O(k„)H3(k)y n ‘ m = 2 M n (m u m 2 ) 


+Vg (( A ‘) i ^)) E E v idu3i) v Ui2,h) 

[\rnij \m 2 J J . _ 


h+ji=k i2+j2=k 


X 


N min(ii ,m\ — 12 ) 

Ni - k Q \ 


m 1 


Z^ 

=0 


A^TVi, 777,1^2) An^mi, /c 0 ,ii, zq) 


x 


At 2 — /co 
m 2 - k 0 


min(j 1 ,m. 2 -j 2 ) 


^2 A U2 (N 2 ,m 2 , j 2 ) Y u (N 2 ,m 2 ,ko,ji,v 2 ) 


122=0 


(70) 

The functions X u and Y\ \ are defined in Eq. (66). Following the same 
procedure as above, the formula for M± 3 is, 


M 13 (m 1 ,m 2 ) = (Ot(k 0 )H3(k)O(k 0 )H(k)) m = 2 (H*(k)) 


2( U\\ rn k+ko,m2-k Q 






(h Z' 2 )}" E E VSM) 

I \m\J \m 2 J I 

h+ji=k i2+j2=k 


X 


x min(i 2 ,mi+ko—ii) 
N i - fc 0 \ 


m i 


Z^ 

!2l=0 


A. ul (Ni,mi + /c 0 ,ii) Xn(iVi,mi,/coZ2,^i) 


x 


N 2 — k 0 
rn 2 - k 0 


min (j 2 ,m 2 —ko—ji) 

Y A U2 (N 2 , m 2 - k 0 ,ji) Y u (N 2 , m 2 , k 0 ,j 2 , u 2 ) 

122=0 


(71) 

Formula for M 22 (mi,m 2 ) is more complicated and we will turn to this now. 
Using the H decomposition H{k ) = J2i+j=k H\{i)H 2 (j) mentioned just 



after the M 11 (mi,m 2 ) formula, the M 2 2 (mi,m 2 ) can be written as, 


M 22 (m 1 ,m 2 ) = (Oi(k 0 )Hi(k)0(ko)Hi(k)) mi ’ m2 

= E E E E Eavg{(O'(Mffi(ii)ft0i)ffiN 

il+jt=k l2+32=k i3+j3=k l4+j4=k 

Xi/ 2 (j 2 )O(A: 0 )i/ 1 fe)i/2(j3)^l(*4)^2(j4)) mi ’ m2 } 


= E E [ Esv s{(OHk 0 )H 1 {i 1 )H 2 (j 1 )H 1 (i 1 )H 2 (j 1 )O(k 0 )H 1 (i 2 )H 2 (j 2 )H 1 (i 2 )H 2 (j 2 )) ,ni ’ m '‘} 

h+ji=k h+j2=k 


Eavg{(Ot(/ Co )// 1 (i 1 ) J ff 2 (ji) J ffi(i 2 )// 2 (j 2 )0(/co)//i(i2) J ff 2 02)^i(h)ff2(ji)> mi ’ m2 } 


+ Eavg{<o t (fc o )ff 1 (i 1 )ff 2 (i 1 )ffi(i2)hr 2 (i2)0(A ; o)^i^i)ff2(ii)hri(f2)ff2O2)> mi ’ m2 }] . 

(72) 

Thus, the ensemble averaged M 22 decomposes into three terms. The first 
term is simple and the next two terms can be decomposed into averages in 
mj and m 2 spaces. Then we have, 

M 22 (m i, m 2 ) = {O^k 0 )O(k 0 )) mum2 {H 2 (k)) mi ' m2 ( H 2 (k)) mi+ko ' m2 - k ° 


+E 2 


o 


Ni 

mi 


N 2 

m 2 


-l 


E E 

h+ji=k i2+32=k 


Y^(i 2 ,j 2 ) 


x X 22:a (Ni,mi, i\,i 2 , k 0 ) Y 22:a (N 2 ,m 2 , j u j 2 , k 0 ) 


+V ( 2 


e> 


iVA /7V 2 
mj \m 2 


E E rg(ii.ii) vg(i 2 ,j 2 ) 

h+ji=k h+j2=k 


x X 2 2 :b(N 1 ,m 1 ,i 1 ,i 2 ,k„) Y 2 2 :b(N 2 ,m 2 ,j 1 ,j 2 ,ko) ■ 

(73) 

Note that the three terms in Eq. (72) correspond directly to the three terms 
in Eq. (73). The second term involves X 22:a and Y 22:a functions. The A" 22;a 
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function is 

X 22 :a(Ni,mi,ii,i 2 ,k 0 ) = ^2 {mi,ax I A a (k 0 ) I rrii + k 0 ,a 2 ) 

ai ,a2 ,a3 ,<24,a5 ,<*6 :a,^i ,1/2 ,Wi/ 2 

x (mi + fc 0 , «4 | ^(fc 0 ) | mi, a 5 ) (mi + fc 0 , a 2 | B” 1 '^ (ii) | m 1 + k 0 , a 3 ) 
x (mi,a 6 | £? I/1AV i(ii) | mi, aq) (mi + k 0 l a 3 \ B V 2 ’ w ^{i 2 ) \ mi + k 0 , a 4 ) 
x (mi,a 5 | B V2 ' tJJu z (i 2 ) | mi,a 6 ) 

( 7 \r 7 \ H ^2 

“EE + /c 0 , mi + fc 0 - ii) A Ul {N 1 , m 1 , m 1 - ii) 

mi / ^=01/2=0 


x A I/2 (iV 1 ,m 1 + k 0 ,mi + k 0 - i 2 ) A U2 (Ni,mi,mi - i 2 )] 1/2 


x 


E 


C 


J'lMa 

/mj+Zeg^ 7 /m^+feg a 3 


c 


/m^6 7 /mj^l 


C 


1^2 i^u<2 

/m^ + fcQ^3 7 /m^+fcQ^4 


s~iko,Oi s~iko,a 

X O „ -7-- w „ -7-- ^ -7-- • 

Jm-^OLS , Jm^OLQ 1 7 /m^+fcQ Q! 2 /m^+/cQ a 4 7 Jm-^OL 5 

(74) 

We have applied Eq. (7) and the Wigner-Eckart theorem to get the second 
form in Eq. (74). Simplification of the six C-G coefficients will finally give a 
compact formula in terns of the X u functions introduced earlier, 





h J 2 

x ^ X u (N 1 ,m 1 ,ko,i 1 ,i' 1 ) y Xu(Ni,m u /c 0 ,i 2 ,^ 2 ) . 

ui=0 U2=0 


(75) 


Similarly, we have 

Y22:a(N 2 ,m 2 ,jl,j2,k 0 ) 


/iV 2 \l 1 / N 2 -k 0 \ 
\koJj \m 2 - k 0 J 


h 32 

x ^ Y u (N 2 , m 2 , k 0 ,ji, zq) ^ y u(N 2 , m 2 , k 0 ,j 2 , v 2 ) ■ 

v\ =0 U2=0 


(76) 
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Finally, the third term in M 22 (rrii,m 2 ) involves the functions A" 22:fe and Y 22:b . 
The expression for function X 22: b is, 

X 22 :b(N 1 ,m 1 ,i 1 , i 2 , ko) = {m i, oq | A a (k 0 ) \ m 1 + k 0 , a 2 ) 


X (mi + k 0 ,a 4 | Al(k 0 ) | mi,a 5 ) (mi + A: 0 ,a 2 | | mi T + aq) 

x (mi, as | B vl ^(ii)B V2 ’ Uv ^(i 2 ) | mi, o;i) . 

(77) 

Now, simplifying Eq. (77) and similarly, Y 22: b, finally we get 

1 /2 

X 22: b(Ni,mi, ii,i 2 , k 0 ) = ^ | : ^)| 

%\ 22 21+22 


x EEEE (mi + fc 0 || [ J B z/1 (ii) J B 1/2 (i 2 )] I/:p || m 1 + k 0 ) 


i/i=0 1/2 =0 i/=0 p 


x (mill [S I/1 (ii)5 V2 (i 2 )] l/:p || mi) (-l)W« 1 +^-/»> 1 . fc o)+Wm 11 /». 1 ,'') 


^ ^v/2721+fco Jmi Jmi+ko Jmi j /fco ^7 ? 

YMN 2 ,m 2 J u h,k„)=(^-_ k ^ {(^ 2 ) d(JV 2 :0} 

jl 22 il+J2 

x EEEE (m 2 II [ir+Oi^C+r || m 2 ) 

i/i =0 1/2=0 i/=0 p 

x (m 2 - fc 0 II [EC 1 (ji)B U 2 (j 2 )} vp II m 2 - fc 0 ) 


X (_!)«(/«*,/m 2 -fc 0 ,fco)+^(/m 2 -fc 0 ,/m a -fc 0 ,*') U (f m2 f m 2 -k 0 fm 2 fm 2 -k 0 \ fk 0 v) . 

(78) 

Combining Eqs. (78), (75), and (76) with (73) will give the formula for 
M 22 (mi,m 2 ). Note that the reduced matrix elements in Eq. (78) are given 
by Eq. (35) along with Eq. (7). 
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6. Asymptotic results for two types of spinless fermions with beta 
and double beta decay type transition operators 


Employing the formulas derived in Section 5 , for some typical values, 
appropriate for atomic nuclei, for N { , A 2 , rn i and m 2 with k = 2 and k 0 
taking values 1 and 2 , numerical results for £(? 77 i,? 77 2 ) and k rs {mi,m 2 ) with 
r + s — 4 are shown in Tables 2 and 3. In the examples shown in the Table, 
N\ and A 2 are sufficiently large but not mi and m 2 for all the examples. For 
the situation with m\ and m 2 sufficiently large (typically larger than say 6), 
the \k rs (mi,m 2 )\ are < 0.3 pointing that the bivariate transition strength 
density is close to a bivariate Gaussian. Also, typically £ ^ 0.6 shows that 
for the systems considered, EGUE is essential (as shown in [ 24 ], £ = 0 for a 
GUE). For further understanding the smallness of k rs (rn 1 , m 2 ), we will derive 
asymptotic formulas first for M rs {mi,m 2 ) and using them, for k rs (mi,m 2 ) 
with r + s — 4 and also for £(7771, m 2 ). 

Let us begin with the asymptotic limit defined by Ni —* 00, N 2 00, 
mi, m 2 fixed with k and ko much smaller than mi and m 2 . Note that in 
the dilute limit (or true asymptotic limit) we also have 7771 —> 00, m 2 —* 00, 
mi/Ni —> 0 and rn 2 /N 2 — > 0 and we will consider this in the later part of 
this section. First, we consider the following functions, 


Ai(A, 777, i, t ) 



d(N : i ) A* (A, m, m 


1/2 


i) A l (N, m + t,m + t — i ) 


x 



fmi fm+ti fm i fti i) 




A 2 (N, 777, i, t) 



d{N : 1 ) A l (N, 777,777 


1/2 


i ) A*(AT, 777 — t,m — t — i) 


x | A\ fN t\ \u( f f f fm _ t . f 

. 777 / V 777 — t / 1 1 


T(N, 777, 7) = A 0 (A, 777, i) , 


A' 


-1 


F(A, 777, i, j) = A*(A, 777,777 — i) A®(A, 777, j) d(N : i) 

V 777 / 


(79) 
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Using Eq. (7) for A, Eq. (14) for d(N : i ) and Eq. (67) for the [/-coefficient, 
Ai and A 2 in the asymptotic limit simplify to 


Ai(N,m,i,t ) = 

a syrn p 

and 

A 2 (N,m,i,t ) = 


N — i + 1 
i 


m'j ^N + 1^ A N — ^.Y — m — i 

m\ /A r \ AN — m — i 


m\ /A[+l\ AN — m\ An — A AN — i + 1 


(80) 


-1 


asyrnp f m\ / N\ Am — i 

i / V i / V t 


Similarly T and F are given by, 


(81) 


, asyrnp f m\ f N 


T(N , m, i 


„,, T . asyrnp f m\ (m — i\ (N\ (N 


(82) 


Using Eqs. (61) and (82), 

£ V^i,j)T(N lt m I ,i)T(N 2 ,m 2t j) . (83) 

Similarly for Mu (mi, m 2 ) given by Eq. (66), in the asymptotic limit only the 
terms with U\ = % in Adi and v 2 = j in Yu will contribute. Then, simplifying 
using Eqs. (80) and (81) we have 

M 11 (m 1 ,m 2 ) aS n? p vg Y. 


i+j=fc 


fco 


fc 0 ; 


T(N 1 ,m 1 ,i)T(N 2 ,m 2 ,j) . 

_____ (84) 

Turning to fourth order moments, firstly for (H 4 (k))" ll ' m2 , using Eqs. (62) 
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and (82), we get 


{H*(k)) mum2 aS -^5 P 2 [ (H 2 (k)) 


i mi,m2 


E E VH(iuhM(%2,h) 

h+ji=k h+j2=k 


x F(JVi, mi ,ti,ij) F(N 2 ,m 2 ,j 1 ,j 2 ) • 

(85) 

This formula follows from the fact that in Eq. (62), in the asymptotic limit 
only terms with U\ = i in X and u 2 = j in Y will contribute. Equation (85) 
also gives M 40 (mi,m 2 ) and M 04 (mi,m 2 ) via Eq. (63). 

The first non-trivial fourth order moment is M 31 (mi,m 2 ) and it is given 
by Eq. (70). As only terms with V\ = i 4 in X n and u 2 = j\ in Y n in Eq. 
(70) will contribute in the asymptotic limit, we have 

M 31 (m u m 2 ) = 2 (H*(k)) mi ™ M 11 (m 1 ,m 2 ) 


E E ^ff(* 1 ,7 1 )^ff(* 2 ,i 2 ) { Ail (ATi, m 4 , i 2 )Ai(Ni,m 1 ,ii, k 0 ) 

ii+ji=k i2+j2=k 


+ A J1 (N 2 , m 2 ,j 2 )A 2 (N 2 , m 2 ,ji, /c 0 )} . 

( 86 ) 

Substituting Eqs. (80) and (81) for A\ and A 2 respectively and Eq. (7) for 
A will give, 

M 31 (m 1 ,m 2 ) aS ^? P 2 (H\k)) mm M u (m u m 2 ) 


+VS E E V N(‘iJi)V3(i 2 J 2 ) 

il+jl=k i2+j2=k 



f m 2 - jA 

V k 0 J 


(87) 


x F(Ni,mi,i 1 ,i 2 ) F(N 2 ,m 2 ,j 1 ,j 2 ) . 

Moment M 13 (m 4 ,m 2 ) is given by Eq. (71) and in the asymptotic limit only 
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the terms with zq = i 2 in X n and u 2 = j 2 in Y n will survive. Then, 


M 13 (m 1 ,m 2 ) = 2 ( H 2 (k )) 


m\+ko,m2—ko 


M u (m 1 ,m 2 ) 


+V£ ^H('ii,ji)V£(i 2 ,j 2 ) A l2 (N 1 ,m 1 + k 0 ,ii)A 1 (N 1 ,m 1 ,i 2 ,k 0 ) 

ii+ji=k i2+j2=k 

x A J2 (N 2 , m 2 - k 0 , ji)A 2 (N 2 , m 2 ,j 2 , k 0 ) . 

( 88 ) 

Substituting Eqs. (80) and (81) for Ai and A 2 respectively and Eq. (7) for 
A will give, 


Mi 3 (mi,m 2 ) 


asyrnp 


2 {H 2 (k)) 


mi+ko,m2—ko 


M n (mi,m 2 ) 


+vg x E 

h+ji=k h+j2=k 


m 1 + k 0 - i 2 

h 


x 


Ni - mi - k 0 + h~ «2\ f mi\ (NX f N x - mi - i 2 




l‘l J \l2 


ko 


m 2 - k 0 - j 2 \ f N 2 m 2 + k 0 + j 1 - j 2 \ /m 2 \ / N 2 \ fm 2 - j 2 \ 

h )\ h )\h)\32)\ k 0 J' 

m 

Finally, the most complicated moment M 22 (m\, rn 2 ) is given by Eq. (73) and 
it has three terms. The first term M^irrii, m 2 ) is simple, 

M<‘>(mi, m 2 ) = (Ot(*„)C >(* 0 )>’” 1 '”' 3 T[P(kff ! ^ . 

(90) 

The second term M^ 2 [mi,m 2 ) is given by Eq. (73) with functions A^ 22:a and 
Y 22 -.a given by Eqs. (75) and Eq. (76) respectively. In the asymptotic limit, 
only the terms with zq = q and v 2 = i 2 will survive in Eq. (75). Similarly, 
only the terms with zq = j\ and zq = j 2 will survive in Eq. (76). These will 
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give 

M$(m 1 ,m 2 ) = V£ Y Y V H( i i,ji)VH(i 2 ,j 2 ) 

il+jl=k i2+j2=k 


X 



(Ni - fc 0 \ 
V m i / 


-i 

Ai(Ni,mi,ii, k 0 ) k 0 ) 


(91) 


x f N2 \f N2 M A 2 (N 2 ,m 2 J 1 ,ko) A 2 (N 2 ,m 2 ,j 2 ,k 0 ) . 

\k 0 J \m 2 J \m 2 - k 0 J 

Substituting Eqs. (80) and (81) for A\ and A 2 respectively and further 

simplifications using the assumption that rrii » k,k 0 and m 2 » k,k 0 will 

give, 


M^(m u m 2 ) as l^ p vg Y V H( i i,ji)VH(i 2 ,j 2 ) 

il+jl=ki 2 +j2=k 


f m 2 - ji - j 2 \ 

V *b ) 


x(* mi k H T(N 1 ,mi,i 1 ) T(N 1 ,m 1 ,i 2 ) T(N 2 ,m 2 ,ji) T(N 2 ,m 2 ,j 2 ) . 

(92) 

/ o\ 

Lastly, the third term M 22 (rn\,m 2 ) of M 22 {mi,m 2 ) is given by Eq. (73) 
along with Eq. (78). In the asymptotic limit only the terms with U\ — L, 
v 2 = i 2 and v — ii + i 2 will survive in X 22:b given by Eq. (78). Similarly, 
only terms with U\ = j i, is 2 = j 2 and u — j | + j 2 will survive in Y 22:b . Then, 


M 22 \ j rni,m 2 ) 



(N 2 - fc 0 \ 

\m 2 - k 0 J 


x Y H V H(ii,hM(i2,j2) {diNnh + iJdiN.-.h+h)} 1 ' 2 

il+jl=k i2+j2=k 


x (rrq + ko\\ || m, + k 0 ) (m, || [^(ii)5«(i 2 )] <1+ia || m,) 

X |L(/ mi+fco / mi fmi+ko fmi 5 fk 0 H +^)| ^2 || [-B J1 (jd)5 J2 (j 2 )f 1+J2 || m 2 ^ 
X ^m 2 - /c 0 || [B J ' 1 (ji)B J ' 2 (j 2 )]' , ' 1+J ' 2 II Ul 2 - |{/(/m 2 fm 2 -k 0 fm 2 fm 2 -k 0 ! fk 0 j 1 + h) \ ■ 

(93) 
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Now, simplifying the reduced matrix elements as in Section 4 along with the 
[/-coefficients using Eqs. (80) and (81) will give the final result, 

{m ,\, m 2 ) a5 Z^ 1P l/j V H( i ^Jl) V H < y i 2 ,j 2 ) 

h+ji=k i2+j2=k 


X 


fNi - mi 

V h 


f m 2 - ji - j 2 \ 
\ k 0 ) 


F{N U mi, ii,i 2 ) F(N 2 , m 2 ,j 1 ,j 2 ) ■ 

(94) 


Combining Eqs. (90), (92) and (94) we have, 

M 22 ( mi , m 2 ) as ™ p (ot(* 0 )o( / t,,)>”»''’« 

+1/ o E E ^(ii,ji)v|(i 2 ,i 2 ) 

U+Ji=fc Z2+j2=k 



X 


m 2 


- Ji 

k 0 


3 2 


1 T(7Vi,mi,ii) T(JVi,mi,i 2 ) T(N 2 ,m 2 ,j 1 ) T(N 2 ,m 2 ,j 2 ) 


mi - *i - *2 

k 0 

(95) 

The asymptotic formulas given by Eqs. (83), (84), (85), (87), (89) and 
(95) are identical to those obtained using the asymptotic theory of Mon 
and French as derived in detail in [25] for EGOE. This agreement gives 
a good check of the exact formulas derived in Section 5. They also show 
(as the asymptotic results should be valid for any k and k 0 ) that the term 
(—I)/!--)+/(•••) [/(...) in Eqs. (66) and (78) will be |[/(... )|, i.e. the phase of 
the [/-coefficient [see Eq. (67)] cancels with the phase factor (—l)/(---)+/(•••) 
in these equations. Rewriting Eq. (67) in a form that extends easily to boson 
systems, we have for U 2 , 


rn 2 


-3 1 
k 0 


32 J F(N U 


mi, i\i i 2 ) F(N 2 , m 2 ,ji,j 2 ) 



[V(f m , fp, fm, fp ; fm-p O] 2 = 


04;E(E7) (iv- 2 ^ + 1 ) 
(Tll’J (Jv + i) 


(96) 


We will discuss extension of Eq. (96) to boson systems in Section 8. 

All the formulas given above simplify further in the dilute limit (or strict 
asymptotic limit) defined by N\ —> 00 , N 2 —* 00 , mi 00 , rn 2 —> 00 , 
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rn i / N { —> 0 and m 2 /iV 2 —> 0 and k and k 0 fixed. Also, assuming that 
= Vff independent of i and j, the reduced moments and cumulants 
will be independent of Vq and V^. With these, the dilute limit formulas 
for £ and k rs with r + s = 4 are obtained using Eqs. (37), (38), (58), (83), 
(84), (85), (87), (89) and (95). The results are as follows. Firstly, it is 
easy to see that M P0 = M P0 /M 00 = (H p (k)) mi ’ rn2 and M 0P = M 0P /M 00 = 
(-f7 p (&)) mi+fc °’ m2-fc °; P = 2,4. Using Eqs. (83), (84) , (85) and (58) we have 
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asymptotic formulas for bivariate cumulants £, fc 40 and k 04 , 


£(mi,m 2 ) = 


M u (mi,m 2 ) 


M 00 (mi, rn 2 ) M 20 (m 4 , m 2 ) M 02 (m 4 , m 2 ) 


1/2 


asymp 


m 2 

: o 




Jl+jl=k 


1 / 2 ' 


x E T(Ni, m\ + /c 0 , i 2 ) T(N 2 , m 2 — k 0 ,j 2 ) 

%2+j2=k 


\ - rn 2 -j 


iTjik V / 




x > | 12 T(N 2 ,m 2 , j) , 


, , v M 40 (mi,m 2 ) 

%(m 1 ,m 2 ) = -j—- ttt - 3 

M 20 (mi,m 2 ) 

E E F(A^i, mi, *i, i 2 ) F(iV 2 , m 2 , j u j 2 ) 

asymp u+ji^/c i2+j2=k 


1 2 


E T(N 1 ,m 1 ,i)T(N 2 ,m 2 ,j ) 

L»+j=fc 


1 , 


, / v M 04 (mi,m 2 ) 

%K,m 2 ) = -33 -— - 3 

M 02 (mi,m 2 ) 


E E F(Ni, mi + & 0 , *i, * 2 ) F(N 2 , m 2 - k 0 ,j 1 ,j 2 ) 

asymp * 1 +j 1 =fe * 2 +/ 2 =fe 


1 . 


E + k 0 ,i)T(N 2 ,m 2 - k 0 , j) 

\_i+j=k 

(97) 

Note that the functions T and F are given by Eq. (82). Similarly, Eqs. (87) 
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and (89) will give the formulas for k :i i and k\ 3 respectively, 


k 31 (m 1 ,m 2 ) = - 


M 3 i(m 1 ,m 2 ) 


3/2 

M 2 o(m 1 ,m 2 ) I M 02 (mi,m 2 ) 


— 3£(mi, m 2 ) 


asyrnp . I /m 2 

-dmi,m 2 ) + { ( 

K 0 


E T(N 1 ,m 1 ,i 1 )T(N 2 ,m 2 ,j 1 ) 

Lh+ji=k 


3/2 


X 


E ^T^mi + ^ 0 , 22 ) T(N 2 , m 2 - k 0 ,j 2 ) 

Ll2+j2=k 


1/2 ' 


-1 


x E E ( m2 fco Jl ) F (- /Vi > mi > ii > i 2 )- F (- /V 2 ,m 2 ,ji,j 2 ), 

il+jl=k i2+j2=k ' 2 


/ci 3 (mi,m 2 ) = - 


Mi 3 (mi,m 2 ) 


1/2 r 


M 20 (mi, m 2 ) M 02 (mi, m 2 ) 


1 3/2 


3£(mi,m 2 ) 


asyrnp . . I I m 2 

-?( m i, m 2 ) + <; ( 

K 0 


E T(Ni,mi,ii) T(N 2 , m 2 ,ji) 

ii+jl=k 


1/2 


X 


1 3/2' 


E T(Ni,mi + k 0 ,i 2 ) T(N 2 ,m 2 - k 0 ,j 2 ) 

\j2+j2=k 


-1 


x E E Vfc 0 J2 ) T (' /Vl,mi,i2 ) T (- /V2 ’m 2 >i 2 ) 

*l+ll=fc *2+12=fc ' ' 


X 


/iVi 

\*1 


mi + A; 0 - i 2 

*i 


/A^ 2 \ /m 2 - ko — j 2 

V/1 / V h 


(98) 
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Finally, using (95) we have, 


k 22 (m u m 2 ) = T 


M 22 (mi,m 2 ) 


M 20 (mi, m 2 ) M 02 (mi, m 2 ) 


-2£ 2 (m 1 ,m 2 ) 


^. 2e{mum2) + 


m 2 \ \ ^ 


kj .4- 




u+jl=fc 


- 1 




- 1 (99) 

x ^ T(Ni, mi + ko, i 2 ) T(N 2 , m 2 — ko,j 2 ) > 

*2+j2=fc / 


x E E (” 2 j2 ) m.mhionM,™!,!.) 

*l+jl=fe 12+J2=fc ' / 


x T(N 2 ,m 2 ,ji) T(N 2 , m 2 ,j 2 ) + F(Ni,mi,ii, i 2 ) F(N 2 , m 2 ,ji,j 2 )} . 


Numerical results obtained using Eqs. (97) - (99) are shown in Tables 2 for 
ko = 2 and k = 2 and similarly in Table 3 for ko = 1 and k = 2. It is 
seen that in general |AtPQ(mi, m 2 )| < 0.3 implying that the strength densities 
are close to a bivariate Gaussian. However, as seen from Tables 2 and 3, 
it is necessary to add the corrections due to kpQ. Also, expanding kpQ in 
powers of 1/mi and l/rri 2 using Mathematica, it is seen that all the kpQ with 
P + Q = 4 behave as, for k = 2 and k 0 — 2, 


kpQ — 





+ ... . 


( 100 ) 


Therefore, for m i >> 1 and m 2 « m ’^ 2 , the transition strength density 
approaches bivariate Gaussian form in general. It is important to recall that 
the strong dependence on rrii in Eq. (100) is clue to the nature of the oper¬ 
ator O i.e., 0{ko) |mi,m 2 ) = |mi + ko,rn 2 — ko)- Thus, we conclude that 
bivariate Gaussian form with corrections (see Appendix A for the form with 
Edgeworth corrections) due to kpQ, P + Q = 4 will form a good approxi¬ 
mation for transition strength densities generated by beta and double beta 
decay type operators. 
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7. Lower-order moments of transition strength densities: results 
for particle removal operators 

Particle removal (or addition) operators are of great interest in nuclear 
physics. For example one particle (proton or neutron) removal from a target 
nucleus gives information about the single particle levels in the target and 
similarly, two-particle removal gives information about pairing force [42, 43]. 
Let us begin with m spinless fermions in N sp states and a particle removal 
operator O that removes ko number of particles when acting on a m fermion 
state. Then the general form of O is, 

O = A M ■ (101) 

ao 

Here, A ao (k 0 ) is a k 0 particle annihilation operator and a 0 are indices for 
a k 0 particle state. Note that A ao (k 0 ) transforms as {fk 0 } = {1^^°} with 
respect to U(N) and At (fc 0 ) transforms as {fk 0 }- It important to recognize 
that the O matrices will be rectangular matrices connecting m particle states 
to m — k 0 particle states. In the defining space, the matrix will be a 1 x d, 0 
matrix with matrix elements given by V ao . Note that «o takes do values 
and do = (^). We will represent the O matrix in the defining space by 
GUE implying that the defining space matrix elements V ao are zero centered 
independent Gaussian random variables. Also, the V ao are assumed to be 
independent of the Vij(k) variables in Eq. (2) and therefore also independent 
of the W variables in Eq. (4)] with variance satisfying 

VaV} = Vl • ( 102 ) 

In many particle spaces the O matrices will be d\ x c ?2 matrices connecting 
d\ = (^) number of m-particle states to d -2 = ( m N ko ) number of (m — ko)- 
particle states. Fig. 4 gives an example for the H and O matrices in the 
defining space and in the m particle spaces. Using Eqs. (101) and (102), we 
have 

{OW) m = ^ , (OOi) m = V%( N ~ . (103) 

Similarly, Eq. (15) gives the relations, 

{(OWHP) m = (OW) m (HP) m , (CAHP(D) m = (OW) m (HP) m ~ k ° . (104) 
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We will also make use of Eq. (57) given before. Following the procedure used 
in Section 2, it is possible to derive formulas for the lower order bivariate 
moments of the transition strength densities generated by O defined by Eq. 
(101). Just as in Sections 3 and 5, we will consider the bivariate moments 

M P q = (i OmQOH p ) m (105) 

with P + Q = 2 and 4 (the P + Q = 3 moments are zero as we are using 
independent EGUE representations for O and 77 matrices in many particle 
spaces). 

7.1. Exact formulas for the bivariate moments 
Firstly, Eqs. (104) gives, 

M 20 = (i OW) m (77 2 ) m , M 02 = <« OW) m (77 2 ) m ~ fc ° , n ^ 

M 40 = (i OW) m (77 4 ) m , M 0 4 = (i OUD) m (77 4 ) m - fco . 

Now, Eq. (103) along with Eqs. (11) and (13) will give the formulas for 
M 2 o, M 02 , M 40 and M 04 . Formula for the first non-trivial moment Mu = 
(Ch77(977) m is derived by introducing complete set of states between 0 1 and 
77, PI and O and O and 77 in the trace giving, 

M n (m) = (OmOH) m = 

-1 

Y. (m, V\ | | m - ko, v 2 ) (m - k 0 , v 3 \ O \ m, v 4 ) (107) 

V\,V2,V3,V4 



x (m — k 0 , v 2 | 77 | m — k 0 , v 3 ) (m, u 4 | 77 | m, tq) . 

Using Eq. (101) and applying Eq. (102) along with Eqs. (4) - (7) and the 
Wigner-Eckart theorem will give, 

/ AT\ _1 

Mn (m) = V^Vl 


N 

m 


k 0 

k 0 


Y. [A 1 "(TV, m — ko , m — k 0 — k) A U (N, m, m — k)] 


1/2 


v =0 


X 


E 


C 


Jk 0 , 


c 


Jk 0 , 


c 


V.UJv 




Vl,V2,V3,V4-, or,u v 


fmV 1 , /m — fcQ^2 fm—kQ V S 5 fmV4 fm — kft V 2 5 fm—kQ^3 fm'V4 ? fmVl 


(108) 


43 



Simplifying the four C-G coefficients will give finally, 


M u (m) 


V 2 V 2 

v o v h 




Z u (N,m, ko, k, v) 


(^) d{N : v) A U (N, m , m — k) A U (N, m — k 0 , m — k 0 



X 


ko 


fm fm—ko i 


fko v) 


(109) 

Here, \U\ appears as discussed just after Eq. (95) and Eq. (96) gives the 
formula for U 2 . 

Turning to the fourth order moments, we need M 13 , M 31 and M 22 . As 
O t 7 ^ O , here M 43 7 ^ M 34 [similarly M 40 7 ^ M 04 and M 2 o 7 ^ M 0 2 as seen 
from Eq. (106)]. Following the procedure used for deriving the formula for 
M u (m), we have for M 31 (m) 


M 31 (m) = (O'HOH 3 )™ = 2 (H 2 ) m M u (m) 


4 . V 2 V 2 
4 - v 0 v H 


-1 


E 


^1 ,V2 ,^3 ,t>4,t>5 ,^6,ZX2 ,CJt/ 2 

(m, v 1 | 4 ),(fco) | m - fc 0 , ^2) (m - fco, v 3 | 4ct (fc 0 ) | m, f 4 ) 


( 110 ) 


x (m — k 0 ,v 2 | B ui ’ u ^(k) \ m — k 0 ,v 3 ) ( m,v 5 \ B vl ' Uu ^{k) \ m,v 6 ) 


x (m, u 4 | B U2 ' Uu ^{k) | m,v 5 ) (m,v 6 \ B V2 ' Wv *{k) \ m, iq) . 

Now, applying the Wigner-Eckart theorem, using the results in Section 2 and 
simplifying the resulting C-G coefficients will give, 

M 31 (m) = (OWOH 3 ) m = 2 (tf 2 ) m M n (m) 




— 1 / AT 7 \ , min{k^m—k) 

0 A u (N,m,k) Zn(N,m,ko,k,v) . 


m — Ay 


i /=0 


(HI) 

The function 1 is defined in Eq. (109). Following the same procedure as 
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above, the formula for M 13 is, 


M 13 (m) = (OW 3 OH) m = 2 (. H 2 ) m ~ ko M n (m) 


min(k,m—kQ—k) 

+VSVS (iv'dzil) V A •'(N,m-k 0 ,k)Z u (N,m,k 0 ,k,v). 

u=0 

( 112 ) 

Formula for M 22 follows from the formula given in Section 5 for M 2 2 (Vfti, m 2 ) 
by using the m 2 part appropriately. The final result (with p a multiplicity 
label) is 

M 22 (m) = (O^H 2 OH 2 ) m = (tf 2 ) m (H 2 ) m ~ k ° 


+v£ V , 2 



Z u (N,m,k 0 ,k,u) 


+Vo V% 


N 


m 


-1 


N-k 0 
m — k 0 


2k 


vi=0 V 2=0 v=0 



d(N : u) 


x ( m II [ B vl {k)B v \k)] v:p || m) (■m - k 0 || [B ul (k)B u \k)] v '- p \\m- k 0 ) 

p 

x (_i)W», ; kQ v). 

(113) 

The moments M P q can be converted into reduced (scale free) cumulants k P Q 
that gives information about the shape of the bivariate transition strength 
density. For our purpose the first non-trivial cumulants are the fourth order 
cumulants and they are defined by Eq. (38). The k P Q, P + Q = 4 follow from 
Eqs. (11), (13), (103), (106),(109), (111), (112) and (113). Numerical results 
for some typical values of (TV, m, k, ko) are shown in Table 4. These results 
show that in general |/cpg| ^ 0.3 indicating that the bivariate strength density 
will be close to a bivariate Gaussian. For further confirming this result, we 
will derive asymptotic results for k P Q. 
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7.2. Asymptotic formulas for bivariate moments and approach to bivariate 
Gaussian form 

Here we will consider the asymptotic limit defined by N — y oo with 
m, k and k 0 fixed and k, k 0 « m. Note that in the dilute limit (or 
true asymptotic limit) we also have m —> oo and m/N —>■ 0 with k and 
ho fixed. Firstly, from Section 6 it is easy to see that in the asymptotic limit: 
(■) Cr‘C:9 Zn(N,m,k„,k,k ) -4 (») A °(N,m,k) -> 

(?) (?); (iii) O'* A*(w,m, m - k) A‘(JV.m, *) d(N : *) -f (”) ("-*) (?)? 

Starting with it should be clear that in the asymptotic limit only the term 
with u — k in Eq. (109) will survive. Then, applying (i) and (ii) above we 
have 


= 


Mu (m) 


(m—k\ (m\ 

asymp ( k(t ) ( fc J 


1/2 


M 00 (m) M 20 (m)M 02 (m) 


1 1/2 


) (V) 


-M 1/2 ' 


(114) 


Similarly, for k 40 and /co 4 only th e terms with z/ = k in Eq. (13) will survive 
and then applying (ii) and (iii) above will give, 


ho(m) = - 


k 04 (m) = 


Mio(m) 

M 20 (m) 

M 04 (m) 
M 02 (m ) 


(m—k\ 

asymp ( fc ) 

ff) 


-3 


asymp 


/m—ko—k\ 

V A: 7 


(115) 


- 1 . 


For M 31 , the first term in Eq. (Ill) is trivial and in the sum in the second 
term only the v — k term will survive in the asymptotic limit. Now, applying 
(i)-(iii) above will give the result for fc 3 i(m), 


k 3 i(m] 


asymp 


/ m—k\ fm—k\ 
V k ) V ko ) 



- £(m) = Cm) k m (m) 


(116) 


Similarly ki 3 (m) is given by, 


k 13 (m, 


fm—ko—k\ (m—k\ (m\ 

asymp ( k ){ kn ){ k ) 


1/2 


m— ko\3/2 


(m\ (m—ko\ 

\knJ V k ) 


-£(m) =£(m) k 04 (m) 


(117) 
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Finally, in M 2 2 only the third term in Eq. (113) is complicated. This is sim¬ 
plified using its relation, valid in the asymptotic limit, to £(m) as described 
in Section 4. Following this we have for k 2 2 , 


k22(m) 


asyrnp 


(m\ (m—k\2 tm—2k\ fm—k' 

9 \P(mX I 2 4- ^ k ° ' to ' ' ^ 

9 - ^ m _ ko ^ -T ^ (rn-ko'j 


-2 [£(m)] 2 + 


(m—2k\ 

\ k 0 ) 
(m\ (m—k o\ 
\k 0 ) \ k ) 


+ 


fm — k\ 

l k ). 


(118) 


In the dilute limit with m —> 00 and m/N 0 and expanding the binomials 
in Eqs. (114) to (118), it is seen that to order 1/m the cumulants k rs , 
r + s — 4 will be —k 2 /m (independent of ko) and the correlation coefficient 
^(m) —> 1 — (. kko)/2m. Thus, the cumulants will tend to zero giving bivariate 
Gaussian form. However, as £ —>• 1 as m —> 00 , in practice it is necessary to 
add the k rs , r + s = 4 corrections to the bivariate Gaussian. 

Let us add that the results given in this Section extend easily to particle 
addition operators O t V Q Hf (fc 0 ), acting on a m-particle state gener¬ 

ating m + k 0 particle states, by using the results in Section 5 for the mi part 
appropriately. 

In addition to the three fermionic systems considered so far, for example 
in nuclear physics applications in particular it is also important to consider 
explicitly the parity symmetry. For this situation, we need to consider two 
orbits one with +ve and other with —ve parity (for two types of fermions, 
we will have four orbits). By combining the formulation given in [44] for 
embedded ensembles with parity with the formulation in Sections 3 and 5, 
it is possible to derive formulas for the bivariate moments over spaces with 
fixed-m [or fixed-(mi, m 2 )] and parity. Results for all these extensions will 
be discussed elsewhere. 


8. Lower-order moments of transition strength densities: Results 
for boson systems 

For m spinless bosons in N sp states with a general Avbody Hamiltonian, 
we have BEGUE(/c) [‘B’ here stands for bosonic]. The embedding algebra 
for this system is again U(N). As m boson states should be symmetric 
under interchange of any two bosons, the irrep f m = {m}; the totally sym¬ 
metric irrep of U(N). Similarly, f m = For H a BEGUE(/c) and 
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the transition operator an independent BEGUE(/c 0 ), the results of Section 
3 translate to those of BEGUE(/c) by applying the well known N —> —Y 
symmetry, i.e. in the fermion results replace N by —N and then take the 
absolute value of the final result; see [13, 14] for discussion on this property 
and [45] for explicit derivation of the formulas for boson systems for the mo¬ 
ments of the one and two-point functions in eigenvalues without using the 
N —> —N symmetry. Firstly, it is easy to see that the m boson space di¬ 
mension DB(N, m ) = ( iV+m_1 ) follows from the m fermion space dimension 
Dp(N, m) = (^) by replacing N by —N in (^) and taking the absolute value 
giving in general, 


(N + r^j N-y-N 


— r + s — 1^ 


(119) 


More strikingly, as pointed out in [14], 


A u (N,m, k ) 


bosons 


A UN,m,k) 


|'m — ^—N — m + k — 


m — v 
k 


N + m + v 
k 


( 120 ) 

Here we have used the formula given by Eq. (7) for fermion systems. More¬ 
over, for bosons the irreps v for a k -body operator take the values v = 
0 ,1 ,..., k as it is for fermions but for the fact that they correspond to the 
Young tableaux {2u, v N ~ 2 }. Also, the N —> — N symmetry and Eq. (14) will 
give 


d B (N : v) 


rry-rr.-y 


( 121 ) 


Using Eqs. (119), (120) and (121), it is possible to write the formulas for 
M rs (m), r + s = 2,4, that correspond to the Eqs. (11), (13), (23), (29) 
and (36). As an example, the bivariate correlation coefficient for the system 
considered in Section 3 but for bosons is given by, 


= 


M n (m) 


min(t,m—k ) 

A b(N, m , k ) A^(A^, m, m — t) ds(N : v) 

i /=0 


y/M 2 o(m) M 0 2 (m) 


N + m — 1 
m 


A ° b (N, to, t) A ° b (N, to, k) 


( 122 ) 
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Again, it can be verified that Eq. (122) has correctly the k t symmetry. 
Let us add that for boson systems, asymptotic results correspond to the 
dense limit defined by m —> oo, N —> oo, rn/N —> oo and k, k 0 fixed. 

Two species boson systems are also important in quantum physics [13, 46] 
and for these systems a situation similar to the one we have considered in 
Section 5 is possible. Again, it is possible to apply the N —N symmetry 
but now for both Ni and N 2 using Eq. (119). Then, ( ,Vl+r ) will change 
to ( JVl_r s +s_1 ) and changes to ( iV2_t + u “ 1 ). Similarly A v (Ni,mi,k), 

d(N\ : u), A u (N 2 ,m 2 ,k) and d(N 2 : u) will change to the bosonic A b and 
ds via Eqs. (120) and (121). Also, for example for /e-particle boson creation 
and annihilation operators A^ a (k) and A a (k), it is easy to prove that 

_|_ jyi _ 

r ) • (123) 

m — k ) 


This follows from the relation (52 a A^(k)A a (k)y" = (”'), the Wigner-Eckart 
theorem and the sum rules for C-G coefficients. More importantly, Eq. (123) 
follows also from the fermionic Eq. (57) by applying N —> —N symmetry. 

Going further, as seen from Sections 3, 5 and 7, the bivariate moments 
contain not only the functions A"(./V, m, k ) and d(N : v) but also U (. N ) Racah 
(or U—) coefficients [(see Eqs. (36) and (66)]. In principle the N —» — N law 
applies also to the Racah coefficients by translating the irreps appropriately. 
For example, the fermionic [/-coefficient given by Eq. (96) can be translated, 
using the N —> —N law, for boson systems and the result is 


/fif ( fm■ fp-> fmi fp i fm—pi n) 


N + v - 2 
v 


m — v\ f N + m + v — 1 


p — v 


m — p 


N + to — 1 
P 


N + m — p — 1 

TO — p 


(N + 2u-l) 


(N — 1) 

(124) 


Note that in Eq. (124), f r = {r}, f s = {s^ 1 } and v = {2 u ,u N ~ 2 }. It is 
easy to verify Eq. (124) for v = 0. Using Eq. (124), it is possible to deal 
with boson systems that are similar to the fermionic systems considered in 
Sections 5 and 7. For example, Eqs. (124) and (109) will give the formula for 
the bivariate correlation coefficient for particle removal operator for bosons. 
Full details of the formulas for the bivariate moments for boson systems and 
their asymptotic structure will be presented elsewhere. 
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9. Conclusions and future outlook 

Embedded random matrix ensembles were used with success in the past 
to understand the form of the eigenvalue density of finite quantum systems 
[16] and now there are good applications of the Gaussian form found for 
these systems [19, 47, 48]. Similarly, these ensembles are shown to provide 
the basis for the theory for expectation values of operators [19, 49]. In the 
present paper the focus is turned to transition strengths. 

Employing embedded Gaussian unitary ensembles of random matrices, for 
the first time we have derived exact group theoretical formulas for the sec¬ 
ond and fourth order bivariate moments of the transition strength densities. 
Explicit results for a spinless many fermion system with /c-body Hamiltonian 
and the transition operator a A;o-body operator are presented in Sections 3. 
Similarly, results for a system with two types of spinless fermions and beta 
decay (and double beta decay) type operators are presented in Section 5. 
In addition, results for a particle removal operator are presented in Section 
7.1. The corresponding asymptotic results presented in Sections 4, 6 and 7.2 
respectively and the numerical results from the exact formulas for the bivari¬ 
ate fourth order cumulants as presented in Tables 1-4 clearly show that the 
smoothed transition densities can be very well approximated by Edgeworth 
corrected bivariate Gaussian. Also, values of the bivariate correlation coeffi¬ 
cient (£) shown in Tables 1-4 clearly confirm that the strength distribution 
will be a narrow distribution (unlike for a GOE). Briefly discussed in Section 
8 are some results for boson systems. The formulation and results given in 
this paper hold for a general A;-body Hamiltonian (similarly for the body rank 
of the transition operators considered). However, it should be emphasized 
that two-body Hamiltonians, i.e. EGUE(2)s , are in general more relevant for 
systems such as nuclei and atoms. Clearly, the work presented in the paper 
represents major progress in random matrix theory for smoothed transition 
strengths after the paper by French et al in 1988 [24], 

One gap in the present results is that there is not yet a formula (or a 
good procedure) available for the U(N) Racah coefficients of the type, 

U(fm , Vl, fm, V2 ■ fm, ») 

where f m = {l m } for fermions (and {m} for bosons) and /i = {2 M , l^” 2 ^} 
(for bosons fi = {2 11 , /j n ~ 2 }). These coefficients are needed for the M 22 
moment. In future, it is also important to consider fermions systems with 
spin as these are of direct interest in mesoscopic systems [50]. Embedding 
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algebra for the ensembles with spin for the fermions is U (hi) ® SU (2) with 
SU(2) generating spin and hi = N/2 [51]. It is also possible to consider 
ensembles with more general 17(G) <g>SU(r) embedding [52]; the ensembles 
with r — 1, 2 and 4 are important for fermionic systems and r — 1, 2 and 
3 for bosonic systems (for example r — 3 is appropriate for spinor BEC). 
The ensembles with U(Q) <g) SU(r) embedding [13, 51, 52] are studied so 
far only for one and two-point functions in eigenvalues. Even for these, the 
group theoretical results are incomplete [51, 52] as the Racah coefficients 
needed, for example for the fourth moment of the eigenvalue density, are not 
yet available. Another important extension is to consider ensembles with a 
mean-field one-body term which is more realistic for systems such as atomic 
nuclei and atoms. 

In summary, in this paper a detailed analytical study of transition strengths 
has been carried out, going beyond the results in [24], using embedded ran¬ 
dom matrix ensembles and established clearly that the form of the transition 
strength densities for isolated finite interacting fermion systems will be gener- 
ically a bivariate Gaussian with fourth order cumulant corrections. Fig. 5 
shows a bivariate Gaussian with and without Edgeworth corrections. We 
expect that in the near future the bivariate Gaussian form will be used in 
practical calculations of transition strengths (see [21] for an attempt in the 
past) just as it is being done at present in a systematic manner by the Michi¬ 
gan group [47, 48] for nuclear level densities with good success using the 
Gaussian form given by the embedded ensembles. 
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APPENDIX A 

Given a bivariate distributions p(x,y), its integral over y[x] the marginal 
density p 2 (y)\pi(x)\. The centroids and variances of these are the marginal 
centroids and variances and say that they are (ei, erf) for pi{x) and 
similarly (e 2 , erf) for p 2 (y)- Note that the M 2 o and M 02 in Section 4 are erf 
and erf respectively. Bivariate Gaussian in terms of the standardized 
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variables x — (x — e\)/<j\ and y = (y — e 2 )/cr 2 is given by, 


Vg ( x , y ) 


1 


exp 


x 2 — 2 £xy + y' 2 \ 

2(1-a j ■ 


(A-l) 


Here £ is the correlation coefficient. Thus, a bivariate Gaussian is defined 
by the five variables (ei, e 2 , 04 , cr 2 , £). Now, the Edgeworth (ED) corrected 
bivariate Gaussian including k rs up to r + s = 4 is given by 

Vbiv-ED{x,y) = ll+(^He 3 o{x,y) + -^-He2i(x,y) 


wp He 12 {x,y) + ^He 03 (x,y) 
2 0 


+ (J^^He4 0 (x,y) + ^f-Hesi(x,y) 


-^-He 22 (x,y) + ^11 He 13 (x,y) + ^ He 04 (x,y) 


k 3 o zj t -\ , k 30 k 2 i 

—He 60 (x,y) H-— He 51 {x, y) 


(A-2) 


C 21 ^30^12 


12 


He4 2 (x,y) 


k 3 ok 03 

36 


k\ 2 k 2 \ 
4 


He 33 (x,y) 


k\ 2 k 21 k 03 

8 + 12 


He 2i {x,y) 


ki 2 k 33 

12 


He 15 (x,y) 


+ Y^He 06 (x,y)yj | yg{x,y) . 

Note that for EGUE and BEGUE, by definition the k rs = 0 for r + s = 3. 
However in practical applications, these will be non-zero though expected 
to be small in magnitude. The bivariate Hermite polynomials He mim2 (x,y) 
in Eq. (A-2) satisfy the recursion relation, 

(1 £ )He m i+i,m 2 (xj y) (x (x, y) 


- rn 1 He mi _ 1)m2 (x, y) + m 2 ^He 

mi ,7712 — 1 (d y ) ■ 

(A-3) 
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This can be solved using He 00 (x,y) = 1, He w {x,y) — (x — £y) / (1 — £ 2 ) 
and He 0 i(x,y) = (y - £x) / (1 - £ 2 ). 
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Table 1: Bivariate correlation coefficient (£) and fourth order bivariate cumulants k rs with r > s and r + s = 4 for the system 
considered in Section 5 for double beta decay type transition operators with ko = 2. Results are given for various values of 
number of sp state Ni and N 2 and number of fermions m\ and m 2 in these sp states respectively. The Hamiltonian body 
rank (fc) is chosen to be k = 2. Results are obtained using the formulas given in Section 5. Note that for the M 22 that is 
needed for k 2 2 we have used Eq. (73) with the third term replaced by the corresponding asymptotic formula given by Eq. (94) 
as the [/-coefficient needed for the finite- (N, m) formula is not available. The finite N and m results are compared with the 
asymptotic limit results (these are given in the brackets). See text for further details. 


N x 

rrii 

n 2 

rn 2 

£ 

&40 

k()4 

^31 

ki3 

^22 

20 

8 

20 

8 

0.66(0.76) 

—0.34(—0.24) 

—0.35(—0.24) 

—0.22(—0.18) 

—0.23(—0.19) 

—0.01(—0.18) 

20 

10 

20 

10 

0.68(0.81) 

—0.31(—0.19) 

—0.32(—0.20) 

—0.21(—0.16) 

—0.22(—0.16) 

0.05(—0.16) 

32 

10 

32 

10 

0.74(0.81) 

—0.26(—0.19) 

—0.26(—0.20) 

—0.19(—0.16) 

—0.19(—0.16) 

—0.04(—0.16) 

32 

12 

32 

12 

0.77(0.84) 

—0.23(—0.16) 

—0.24(—0.16) 

—0.18(—0.14) 

—0.18(—0.14) 

—0.01(—0.14) 

32 

16 

32 

16 

0.78(0.88) 

— 0 . 21 (— 0 . 12 ) 

— 0 . 22 (— 0 . 12 ) 

—0.17(—0.11) 

—0.17(—0.11) 

0.06(—0.11) 

32 

10 

44 

8 

0.73(0.78) 

—0.27(—0.22) 

—0.29(—0.23) 

—0.20(—0.17) 

—0.21(—0.18) 

—0.07(—0.18) 

32 

10 

44 

15 

0.79(0.85) 

—0.21(—0.16) 

—0.22(—0.16) 

—0.17(—0.13) 

—0.17(—0.13) 

—0.03(—0.13) 

32 

10 

44 

20 

0.80(0.87) 

—0.20(—0.13) 

—0.20(—0.13) 

—0.16(—0.11) 

—0.16(—0.11) 

O.Ol(-O.ll) 

32 

12 

44 

20 

0.81(0.88) 

—0.19(—0.12) 

—0.19(—0.12) 

—0.16(—0.11) 

—0.16(—0.11) 

0 . 02 (— 0 . 11 ) 

32 

16 

44 

20 

0.81(0.89) 

—0.19(—0.11) 

—0.19(—0.11) 

1 

O 

Oi 

T 

0 

0 

—0.15(—0.10) 

O 

O 

T 

0 

0 

44 

10 

58 

20 

0.83(0.87) 

—0.18(—0.13) 

—0.17(—0.13) 

—0.15(—0.11) 

—0.14(—0.11) 

—0.03(—0.11) 

44 

15 

58 

20 

0.84(0.89) 

—0.16(—0.11) 

—0.16(—0.11) 

—0.13(—0.10) 

1 

O 

T 

0 

0 

-O.Ol(-O.lO) 

44 

20 

58 

20 

0.85(0.90) 

—0.15(—0.10) 

—0.16(—0.10) 

—0.13(—0.09) 

—0.13(—0.09) 

0.01(—0.09) 

44 

10 

58 

24 

0.83(0.88) 

—0.17(—0.12) 

—0.16(—0.12) 

1 

O 

hP- 

T 

0 

0 

—0.14(—0.10) 

-O.Ol(-O.lO) 

44 

10 

58 

28 

0.84(0.90) 

—0.16(—0.10) 

—0.16(—0.10) 

—0.14(—0.09) 

—0.13(—0.09) 

0.01(—0.09) 



Table 2: Bivariate correlation coefficient (£) and fourth order bivariate cumulants k rs with r > s and r + s = 4 for the system 
considered in Section 5 for beta decay type transition operators with fco = 1- Results are given for various values of number 
of sp state N\ and N 2 and number of fermions toi and m 2 in these sp states respectively. The Hamiltonian body rank (k) 
is chosen to be k = 2. Results are obtained using the formulas given in Section 5. Note that for the M 22 that is needed 
for k 2 2 we have used Eq. (73) with the third term replaced by the corresponding asymptotic formula given by Eq. (94) as 
the [/-coefficient needed for the finite-(iV, to) formula is not available. The finite N and to results are compared with the 
asymptotic limit results (these are given in the brackets). See text for further details. 



mi 

n 2 

m 2 

■e 

k±o 

ko4 

^31 

ki3 

k 2 2 

20 

8 

30 

8 

0.83(0.88) 

—0.32(—0.24) 

—0.33(—0.25) 

—0.26(—0.21) 

—0.27(—0.22) 

—0.12(—0.22) 

20 

8 

30 

10 

0.84(0.89) 

—0.30(—0.22) 

—0.30(—0.22) 

—0.25(—0.19) 

—0.25(—0.20) 

—0.10(—0.20) 

20 

10 

30 

10 

0.84(0.90) 

—0.29(—0.20) 

—0.30(—0.20) 

—0.24(—0.18) 

—0.25(—0.18) 

—0.07(—0.18) 

20 

10 

30 

15 

0.86(0.92) 

—0.26(—0.16) 

—0.26(—0.16) 

—0.22(—0.14) 

—0.23(—0.15) 

—0.03(—0.15) 

36 

8 

36 

8 

0.85(0.88) 

—0.29(—0.24) 

—0.29(—0.24) 

—0.24(—0.21) 

—0.25(—0.21) 

—0.15(—0.21) 

36 

12 

36 

12 

0.88(0.92) 

—0.22(—0.16) 

—0.22(—0.16) 

—0.20(—0.15) 

—0.20(—0.15) 

—0.08(—0.15) 

36 

16 

36 

16 

0.90(0.94) 

—0.20(—0.12) 

—0.20(—0.12) 

—0.18(—0.12) 

—0.18(—0.12) 

—0.03(—0.12) 

36 

18 

36 

18 

0.90(0.94) 

—0.19(—0.11) 

—0.19(—0.11) 

—0.17(—0.10) 

—0.17(—0.10) 

-O.Ol(-O.lO) 

36 

8 

36 

10 

0.86(0.89) 

—0.27(—0.22) 

—0.27(—0.21) 

—0.23(—0.19) 

—0.23(—0.19) 

—0.13(—0.19) 

36 

8 

36 

12 

0.87(0.90) 

—0.25(—0.20) 

—0.25(—0.19) 

—0.22(—0.18) 

—0.22(—0.17) 

—0.11(—0.17) 

36 

8 

36 

16 

0.88(0.92) 

—0.23(—0.16) 

—0.23(—0.16) 

—0.21(—0.15) 

—0.20(—0.15) 

—0.08(—0.15) 

44 

8 

44 

8 

0.85(0.88) 

—0.28(—0.24) 

—0.28(—0.24) 

—0.24(—0.21) 

—0.24(—0.21) 

—0.16(—0.21) 

44 

12 

44 

12 

0.89(0.92) 

—0.21(—0.16) 

—0.21(—0.16) 

—0.19(—0.15) 

—0.19(—0.15) 

—0.09(—0.15) 

44 

20 

44 

20 

0.91(0.95) 

—0.16(—0.10) 

—0.16(—0.10) 

—0.15(—0.09) 

—0.15(—0.09) 

—0.02(—0.09) 



Table 3: Bivariate correlation coefficient (£) and fourth order bivariate cumulants k rs 
(= k sr ) with r > s and r + s = 4 for various values of number of sp state ( N ), number 
of fermions (m), Hamiltonian body rank (k) and the body rank (ko) of the transition 
operator. Results are obtained using the formulas given in Section 3. Note that for the 
M 22 that is needed for £ 22 , we have used Eq. (36) with the third term replaced by the 
corresponding asymptotic formula given by Eq. (48) as the 17-coefficient needed for the 
finite-(IV, in) formula is not available. The finite N and m results are compared with the 
asymptotic limit results (these are given in the brackets). See text for further details. 
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k 

k 0 
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/t40 

^31 

^22 

20 

10 

2 

1 

0 .68(0.8) 

—0.54(—0.38) 

—0.36(—0.3) 

—0.09(—0.27) 

30 




0.73(0.8) 

—0.48(—0.38) 

—0.35(—0.3) 

—0.16(—0.27) 

50 




0.76(0.8) 

—0.43(—0.38) 

—0.33(—0.3) 

—0.21(—0.27) 

80 




0.78(0.8) 

—0.41(—0.38) 

—0.32(—0.3) 

—0.23(—0.27) 

40 

8 

2 

1 

0.71(0.75) 

—0.52(—0.46) 

—0.37(—0.35) 

—0.23(—0.3) 


12 



0.78(0.83) 

—0.4(—0.32) 

—0.31(—0.27) 

—0.15(—0.24) 


15 



0.81(0.87) 

—0.36(—0.26) 

—0.29(—0.22) 

—0.1(—0.21) 


20 



0.82(0.9) 

—0.32(—0.2) 

—0.27(—0.18) 

—0.03(—0.17) 

60 

8 

2 

1 

0.72(0.75) 

—0.5(—0.46) 

—0.36(—0.35) 

—0.26(—0.3) 


15 



0.83(0.87) 

—0.32(—0.26) 

—0.26(—0.22) 

—0.15(—0.21) 


20 



0.86(0.9) 

—0.27(—0.2) 

—0.23(—0.18) 

—0.09(—0.17) 


30 



0.88(0.93) 

—0.23(—0.13) 

—0.2(—0.12) 

—0.02(—0.12) 

24 

8 

2 

1 

0.67(0.75) 

—0.56(—0.46) 

—0.38(—0.35) 

—0.18(—0.3) 



2 

2 

0.44(0.54) 

—0.56(—0.46) 

—0.24(—0.25) 

—0.07(—0.17) 



3 

1 

0.54(0.63) 

—0.88(—0.82) 

—0.47(—0.51) 

—0.24(—0.35) 



3 

2 

0.27(0.36) 

—0.88(—0.82) 

—0.23(—0.29) 

—0.06(—0.12) 



4 

1 

0.41(0.5) 

—0.99(—0.99) 

—0.41(—0.49) 

—0.16(—0.25) 

40 

12 

1 

1 

0.89(0.92) 

—0.11(—0.08) 

—0.1(—0.08) 

—0.02(—0.08) 



2 

1 

0.78(0.83) 

—0.4(—0.32) 

—0.31(—0.27) 

—0.15(—0.24) 



3 

1 

0.68(0.75) 

—0.71(—0.62) 

—0.48(—0.46) 

—0.27(—0.37) 



3 

2 

0.45(0.55) 

—0.71(—0.62) 

—0.32(—0.34) 

—0.12(—0.21) 



4 

1 

0.59(0.67) 

—0.91(—0.86) 

—0.53(—0.57) 

—0.3(—0.4) 



4 

2 

0.33(0.42) 

—0.91(—0.86) 

—0.3(—0.36) 

—0.1(—0.17) 



5 

1 

0.5(0.58) 

—0.98(—0.97) 

—0.49(—0.57) 

—0.24(—0.34) 



5 

2 

0.24(0.32) 

—0.98(—0.97) 

—0.23(—0.31) 

—0.05(—0.1) 


59 






Table 4: Bivariate correlation coefficient (£) and fourth order bivariate cumulants k rs 
with r + s = 4 for various values of number of sp states (N), number of fermions ( m ), 
Hamiltonian body rank ( k ) and the rank (fco) of the particle removal transition operator. 
Results are obtained using the formulas given in Section 7.1. Note that for the M 22 that is 
needed for ^ 22 , we have used Eq. (113) with the third term replaced by the corresponding 
asymptotic formula given by Eq. (118) as a formula for the reduced matrix elements in 
Eq. (113) is not available. 


N 

m 

k 

ko 

£ 

^40 

ko4 

k:i 1 

ki3 

^22 

20 

10 

2 

1 

0.82 

-0.54 

-0.55 

-0.44 

-0.45 

-0.21 

30 

10 

2 

1 

0.85 

-0.48 

-0.50 

-0.41 

-0.43 

-0.26 

60 

10 

2 

1 

0.88 

-0.42 

-0.46 

-0.37 

-0.40 

-0.30 

80 

10 

2 

1 

0.88 

-0.41 

-0.45 

-0.36 

-0.39 

-0.31 

50 

12 

2 

1 

0.89 

-0.38 

-0.40 

-.034 

-0.36 

-0.25 


15 

2 

1 

0.91 

-0.33 

-0.35 

-0.30 

-0.31 

-0.19 


20 

2 

1 

0.92 

-0.29 

-0.29 

-0.26 

-0.27 

-0.13 


25 

2 

1 

0.92 

-0.27 

-0.27 

-0.25 

-0.25 

-0.08 

24 

8 

2 

1 

0.82 

-0.56 

-0.61 

-0.46 

-0.49 

-0.31 



2 

2 

0.66 

-0.56 

-0.67 

-0.37 

-0.43 

-0.22 

40 

15 

2 

1 

0.90 

-0.36 

-0.37 

-0.32 

-0.33 

-0.18 



2 

2 

0.80 

-0.36 

-0.38 

-0.29 

-0.31 

-0.12 

60 

20 

2 

1 

0.93 

-0.27 

-0.27 

-0.25 

-0.25 

-0.14 



3 

1 

0.89 

-0.51 

-0.53 

-0.46 

-0.47 

-0.30 



3 

2 

0.79 

-0.51 

-0.54 

-0.40 

-0.43 

-0.22 


60 
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Figure 1: Schematic figure showing transition strengths, (a) Transition strengths for tran¬ 
sitions induced by an operator 0\, from levels with energies (eigenvalues of the Hamilto¬ 
nian) Ei of a system ‘a’ to levels with energies Ef of the same system ‘a’. The strengths 
| (a,Ef | Oi | a,Ei) | 2 are proportional to the widths of the lines in the figure, (b) Transi¬ 
tion strengths for transitions induced by an operator O 2 , from one particular level with 
energy (eigenvalue of the Hamiltonian) Ei of a system ‘a’ to levels with energies Ef of 
another system ‘b\ The strengths \{b,Ef \ 0\ \ a,Ei) | 2 are proportional to the widths of 
the lines in the figure. In general, transitions from several levels of the system ‘a’ to the 
levels of system ‘6’ are possible as shown in the inset figure. 
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U(Q) Young Tableaux {/}, /i > f 2 > • • • > fn > 0 



(a) general irrep 


If} = {8} 


(b) symmetric irrep 


{f} = {l 7 } 


(c) anti-symmetric irrep 


{f} = {5,3 2 ,2, l 2 } {f} = {6,4,3, l 2 } 



(d) conjugate irrep {f}, given the irrep {f} 



Figure 2: Young tableaux representation of the irreps of U(Sl). Shown are examples of: (a) 
a general irrep {/}; (b) symmetric irrep {m}; (c) antisymmetric irrep {l m }; (d) conjugate 
irrep {/} that corresponds to a given {/}; (e) irrep {/} that corresponds to a given {/}. 
Note the importance of in defining {/}. 
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N x = 20, N 2 = 16, k = k 0 = 2 


r/(2.0) 17(2,0) 
n 21 n 22 

0 0 ••• 0 

0 0 ••• 0 

0 0 

0 0 0 


0 0 ••• 0 

0 0 0 

0 

sji’ 1 * H$ l] - 

ff (M) ff (M) 
n 21 n 22 

0 0 0 

0 0 0 



0 0 0 

0 

0 

c 2) ••• 

w (0,2) 77(0,2) 
n 21 n 22 


} ( 2 , 0 ), 190 


( 0 , 2 ), 120 

(a) H matrix in the defining 
space 


ds 


120 


On 0 12 Oi 3 
Oil Oil O 23 
O31 O32 O33 


190 


[ On 

Ol2 

0 :3 .1 


O21 

O22 

023 . 











rr(6,8) 
-“11 

o-(6,8) 
“12 

rr(6,8) 
“13 

rr(6,8) 
“21 

rr(6,8) 
“22 

rr(6,8) 
“23 

tt{ 6,8) 
“31 

tt( 6,8) 
“32 

tt( 6,8) 
“33 


(b) O matrix in the 
defining space 


„(8,6) „(8,6) rr(8,6) 

“ll “12 “13 

if(M) 7,(8, 6) „(8,6) 
“21 “22 “23 

„(8,e) w (8,6) „(8,6) 

“31 “32 “33 


> cU 


(c) H matrix in the 

= (6,8) space 
dA 


Oli Oli o( 3 
oh O | 2 o | 3 
oil ol o\ 3 


(d) O matrix connecting (6,8) to 
(8, 6) space 


(e) H matrix in the 

(mi,m 2 ) = (8,6) space 


> ds 


(f) O ' matrix connecting 

(8,6) to (6,8) space 


Figure 3: Schematic diagram showing the matrix representations of the Hamiltonian ( H) 
and the transition operator ( O ) in the respective defining spaces and in the many particle 
spaces for the system considered in Section 5 with N\ = 20, N 2 = 16, k = 2 and kg = 2. (a) 
H matrix in the defining space. Note that the matrix is in block diagonal form with 3 blocks 
and they correspond to (mi,m 2 ) = (2,0), (1,1) and (0,2) as shown in the figure. The 
matrix dimensions of each block are also shown in the figure, (b) O matrix in the defining 
space. Note that O is a rectangular matrix connecting (mi,m 2 ), = (0,2) to (mi,m, 2 )/ = 
(2, 0). The number of rows and columns in the matrix are shown in the figure, (c) H matrix 
Hi in the initial (mi, m 2 ) = (6,8) space with matrix dimension = 498841200. (d) O 
matrix connecting states in (mi, 1712)1 = (6,8) space with the states in (rni,m 2 )f = (8,6) 
space with matrix elements {(mi,m 2 )f = (8,6)/? | O | (mi,m 2 )i = (6,8),a). Note that 
the matrix is a ds x d^ rectangular matrix with dA given in (c) and ds = 1008767760. (d) 
H matrix Hf in the final (8,6) space, (e) same as (d) but for the Ch matrix. It is useful 
to note that the bivariate moments Mpq of the transition strength density are given by 
(d^) _1 TV [O* 1 (Hf)Q0(Hi) p ] where TV stands for the matrix trace. 
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N = 20, k = 2, k 0 = 1 


190 


r 

11 a 

rr( 2 ) 
11 \2 

/t( 2 ) . 

^13 


rr( 2 ) 

^21 

rj( 2 ) 

n 22 

/f( 2 ) . 

n 23 


rr( 2 ) 

“31 

o-( 2 ) 

32 

/t( 2 ) . 

11 33 







(a) H matrix in the defining 
space 


20 


On 0 12 . 

(b) O matrix in the 
defining space 


d>A 


rr(10) rr(10) rr(10) 

W 11 -“12 -“13 

|j(10) ej-(IO) ry(10) 

JJ 99 




(c) H matrix in the m = 10 
space 


d B 


On O12 O13 
O21 O22 O23 
O31 O32 O33 


> d A 


(d) O matrix connecting 
m = 10 to m = 9 space 


d,B 


HiV Hif 

uC 9 ) rr( 9 ) H^ 
11 21 n 22 n 22> 

Hfff 


ds 


(e) H matrix in the 
m = 9 space 


d A 


Oil Ot 2 0+3 
oii o+ 2 0+3 


d,B 


(f) O ' matrix connecting m = 9 to 
m = 10 space 


Figure 4: Schematic diagram showing the matrix representations of the Hamiltonian ( H) 
and the transition operator ( O ) in the respective defining spaces and in the many particle 
spaces for the system considered in Section 7 with N = 20, k = 2 and ko = 1. (a) H matrix 
in the defining space, (b) O matrix in the defining space, (c) H matrix Hi in the initial 
to = 10 space with matrix dimension dA = 184756. (d) O matrix connecting states in nrii = 
10 space with the states in rrif = 9 space with matrix elements (to = 9, /3 | O \ m = 10, a). 
Note that the matrix is a ds x dA rectangular matrix with dA given in (c) and ds = 167960. 
(e) H matrix Hf in the final m = 9 space, (f) same as (d) but for the dfi matrix. See Fig. 
3 for further details. 
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( b ) Bivariate Gaussian with ED corrections 


Figure 5: Schematic diagram showing (a) normalized bivariate Gaussian and (b) normal¬ 
ized bivariate Gaussian with Edgeworth (ED) corrections (see Appendix A). The fourth 
order cumulants used are £ = 0.83, = —0.18, fco 4 = —0.17, & 3 i = —0.15, k\% = —0.14, 

and fc 22 = —0.03 corresponding to = (44,10) and (_/V 2 , m 2 ) = (58, 20) in Table 2. 
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